QUADRO 1 
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pára-raios. Esta tensão é dada muito aproximada- 
mente por 


Up= rpi= KW (9) 


No que segue admitimos y = 0,25 o que é válido, 
com boa aproximação, para alguns dos pára-raios 
usados entre nós. 

A solução gráfica das equações 8, apresentada na 
figura 9, supõe u,= O e dá-nos as intensidades de 
crista mais elevadas, imx, das correntes através de 
pára-raios instalados na extremidade de linhas, de 
neutro isolado, a 15 kV, 30 kV e 60 kV. Admite-se 
que os níveis de isolamento das linhas são 200 kV, 
300 kV e 500 kV, respectivamente, para as tensões 
nominais 15 kV, 30 kV e 60 kV. 


Na figura 9 linearizaram-se as características up (1) 
nos troços em que interceptam os segmentos de recta 
Zu; — 1 (Z5, + Rj). Para cada valor de u; (200 kV, 
300 kV e 500 kV) tomaram-se dois valores para Z, 
(200 82 e 400 42) Jedois valores para R;(0e 20 OQ ). 

No quadro | resumem-se os dados que se tomaram 
para elaborar a figura 9 e os resultados que se 
obtiveram. 


Da figura 9 ressalta a importância da impedância 
de onda para a redução da intensidade i através do 
pára-raios. Esta intensidade será tal que o produto 
(Zo + R$) adicionado à tensão up (pouco variável 
com 1) seja igual ao dobro da tensão incidente ui; 
Tem, também, interesse verificar como é reduzida a 
influência da resistência de terra do pára-raios (para 
Ride 0 a 20 ohms) na intensidade que o atravessa. 

Com base no que se expôs até aqui já é possível 
fundamentar uma orientação: “desde que os últimos 
vãos das linhas aéreas estejam eficazmente cobertos 
por cabos de guarda (e a instalação terminal também) 
a Intensidade através dos pára-raios não excederá 5 kA, 
pelo que o nível de protecção nominalmente assegu- 


o+R  HUZ,+R,) 


(i=5 KA) 


] [kV] 


rado por pára-raios de classe 5 KA estará, de facto, 
garantido”. 

Admite-se que os valores de crista das sobretensões 
de manobra não ultrapassam o nível de isolamento 
da linha. Esta hipótese está subjacente à conclusão 
que se tirou. 


** * * 


Diferente será a conclusão se as linhas aéreas esti- 
verem, em toda a sua extensão, desprovidas de cabo 
de guarda e os valores de R, não forem adequados. 

Na figura 10 representa-se o vão, PF P, de uma 
linha aérea, adjacente à instalação terminal. Não há 
cabo de guarda. 

Designa-se por P;o terminal onde estão ligados os 
pára-raios e o transformador protegido e por Pço 
apoio que antecede o terminal. 

Um impacte directo em P; começa por provocar o 
escorvamento do pára-raios, emitindo-se simultanea- 
mente uma onda de tensão u em direcção a P, Se 
admitirmos uq = O as intensidades i (através do pára- 
-raios) e 1, (onda de corrente emitida em direcção a 
P,) são dadas por 


a : Zo 
a Zotrp+R; 


l 


E rp+R; 
CE = TI E = 


em que i, é a intensidade total no ponto do impacte, 


Na maioria dos casos os valore de Z, e R, diferem 
de uma ordem de grandeza e, para altas intensidades 
de descarga, R, e rp também diferem de uma ordem 
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de grandeza. Poderíamos pois concluir que a quase 
totalidade da intensidade ig iria atravessar o pára-raios. 

Todavia se a tensão entre os terminais do pára-raios 
cresce pouco com 1, o mesmo não acontece com a 
tensão entre Pç e a terra, influenciada pelo termo 
linear 1 R,. Esta subida de tensão em Pr(em relação à 
terra) pode provocar novas descargas que estabele- 
cem percursos em paralelo com o pára-raios inicial- 
mente escorvado. Também se deve ponderar que o 
crescimento da tensão u, emitida em direcção a Pc. 
poderá provocar o contornamento da cadeia de iso- 
ladores instalada em P,. 

Contudo, e não obstante os vários percursos de 
descarga, em paralelo, que a subida da tensão vai 
originando (com probabilidades crescendo com Rj), a 
intensidade através do pára-raios pode atingir valores 
elevados [13]. 

A salvaguarda dos pára-raios, quando não há li- 
nhas de guardas reside na fraquissima probabilidade 
da ocorrência de uma descarga directa sobre o pró- 
prio pára-raios ou ao longo do vão imediatamente 
adjacente. Todavia se este vão não estiver eficazmen- 
te protegido não está eliminada a hipótese de uma 
descarga directa atingir Pç Isto é um argumento a 
favor do uso de classe de intensidade superior a SKA 
quando não há linha de guarda. 


9-Janeiro-1981 


a4 


CF. WAGNER, A.R. HILEMAN 


. GW. BROWN, 5. THUNANDER 


FIGURA 10 — Descarga directa sobre o pára-raios 
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APÊNDICE 1 


Quando se dá o contornamento em P, o mode- 
lo é o que se representa na figura 1 1.a. Este mode- 
lo pode simplificar-se 
conforme representa a figura 11.b; i, é uma fonte de 
corrente e não há acoplamento entre os dois ramos 
do condutor activo que divergem de A. 


Então 
= nt; (11) 
Z 
lz 3 h Ky 
£o 
TR, lo—z 


2R, 
— l = 
* Viki+Zs 


1 (12) 


Por cada ramo do condutor activo parte uma onda de 
corrente dada por 


o R, 
2 “MM + 


e a correspondente onda de tensão é dada por 


Zo 


ly; | 
DE RR AR pp 
Mara. 


Tudo se passa como se uma onda de tensão, inci- 
dente, u/ = 1, Rr atingisse o ponto A. A tensão uem A 
está relacionada com u; pelo facto de transmissão K ,5 


£o 


K' E q —— 
0 +. 


FIGURA 11 — Contornamento da cadeia em P,. 


que satisfaz à condição 


u=u;j. K' 


Perante uma onda incidente dada, u; é natural que 
a onda transmitida u cresça com Z,, 
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APLICAÇÃO DA VOLTAMETRIA 
CÍCLICA A CARACTERIZAÇÃO DE 
PROCESSOS DE ELÉCTRODO 


M.A.G. MARTINS e C.A.€C. SEQUEIRA 


Laboratório de Electroquímica, Instituto Superior Técnico, Lisboa 


RESUMO 


Apresenta-se uma compilação dos critérios de diagnoós- 
tico, baseados em estudos teóricos de voltametria ci- 
clica, que permitem caracterizar sistemas de eléctrodo 
desconhecidos. 


ABSTRACT 


ln this paper, it is given a compilation of conve- 
nient diagnostic criteria, based on several theoretical 
treatments of cyclic voltammetry, which enable both 
the qualitative and quantitative characterization of 
unknown electrode processes. 


INTRODUÇÃO 


A caracterização de reacções electroquímicas de 
oxidação e redução pode fazer-se à custa de diversas 
técnicas de entre as quais a voltametria de varrimento 
cíclico, a cronopotenciometria e a cronoamperome- 
tria são as mais vulgarmente utilizadas. Desde 1950 
foram publicados diversos livros (1-4) e artigos de 
sintese (5-7) em que se examinam detalhadamente as- 
pectos teóricos e experimentais dessas técnicas. Nesta 
comunicação, cujo tema é a voltametria cíclica, apresen- 
ta-se uma compilação dos critérios de diagnóstico 
que, sob o nosso ponto de vista, são mais convenien- 
tes para a caracterização qualitativa e quantitativa de 
sistemas de eléctrodo desconhecidos. 


VOLTAMETRIA CÍCLICA 


A voltametria cíclica é uma técnica particularmen- 
te adequada ao estudo de mecanismos de processos 
globais complexos, nos quais intervêm reacções qui- 
micas e/ou electroquímicas, ocorrendo em eléctrodos 
sólidos. 

Foi introduzida em 1938 por Matheson e Nichols 
(8), tendo sido descrita teoricamente com algum 
pormenor, dez anos mais tarde, por Randles (9) e 
Seveik (10). 

Trata-se duma técnica electroquimica em que o po- 
tencial de um eléctrodo indicador ou de trabalho, 
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FIGURA 1 — Relação potencial /tempo, no caso 
da voltametria cíclica 


mantido estacionário numa solução em repouso, que 
contém espécies electroactivas, varia linearmente com 
o tempo. 

Denomina-se voltamograma, o registo da corrente 
em função do potencial aplicado ao eléctrodo, no 
decurso dum intervalo de tempo, em que este poten- 
cial varia periodicamente, de acordo com uma onda 
triangular (com a forma dum triângulo isósceles). 

No caso da voltametria cíclica, este intervalo de 
tempo será, pelo menos, o correspondente a um ciclo, 
envolvendo assim um varrimento anódico e outro ca- 
tódico (Fig. 1). 

As velocidades de variação do potencial aplicado, 
podem oscilar entre alguns mV/h e algumas dezenas 
de V/seg, revestindo-se esta última gama de especial 
interesse em voltametria cíclica no estudo de proces- 
sos de eléctrodo muito rápidos. 

A duração dum ciclo é um dos parâmetros mais 
úteis em voltametria cíclica. Diminuindo esta dura- 
ção, ou seja aumentando a frequência dos ciclos (n.º 
de ciclos / unidade de tempo) de modo conveniente, é 
possível, em certas condições, observar a existência 
de espécies intermediárias com tempos de vida muito 
curtos. 

Para a maior parte dos casos, a informação mais 
importante acerca dum determinado sistema desco- 
nhecido, é extraída dos primeiros cinco ciclos, não 
havendo normalmente vantagem em prosseguir a ex- 
periência para além do 5.º ciclo (2). 

Admitindo que a solução contém inicialmente ape- 
nas a forma oxidada, (Ox), o varrimento de potencial 


Potencial 
aplicado 


E; — Potencial inicial 


E, — Potencial de inversão 


t — Tempo de um ciclo 
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inicia-se a um valor de potencial anódico (mais posi- 
tivo), em relação ao valor do potencial padrão do 
sistema redox (Ox, Red), potencial ao qual, como se 
sabe, não flui corrente através do eléctrodo. 

À medida que o varrimento prossegue e nos apro- 
ximamos do potencial padrão, inicia-se a descarga da 
espécie Ox, com o correspondente aparecimento da 
forma reduzida, (Red). 

A descarga da espécie Ox, é em geral acompanhada 
pelo empobrecimento desta espécie junto do eléctro- 
do (camada pré-electródica), estabelecendo-se um 
gradiente de concentração, entre a concentração de 
Ox na solução, (C), e a concentração de Ox à superfi- 
cie do eléctrodo. À transferência de carga, passa pois 
a associar-se um processo de difusão. 

A valores muito próximos do potencial padrão, a 
corrente cresce bruscamente, até que se atinge um 
máximo, a partir do qual se observa um decréscimo 
da corrente. Este aumento da corrente, devido à reac- 
ção do eléctrodo, é contrabalançado pela deflexão da 
concentração da espécie electroactiva Ox à superfície 
do eléctrodo, proveniente do aumento gradual da es- 
pessura da camada de difusão. A deflexão de Ox por 
conseguinte, origina a queda da corrente. 

Quando se inverte o sentido do varrimento de po- 
tencial, inicia-se o processo de reoxidação da espécie 
reduzida Red, de que resulta um pico anódico. 

A curva corrente potencial (voltamograma cíclico), 
para um processo de redução reversível simples, 
apresenta-se na Fig. 2, com os parâmetros mensurá- 
veis obtidos. 

A corrente de pico (ip) para um processo reversível 
simples, controlado por difusão, é dada pela equação 
de Randles-Sevcik (1): 


p= 0,452 (nF)32(RT) 12 AC (Dy) 12 


em que 


corrente farádica (mA); 

n = número de electrões envolvidos no processo 
de eléctrodo; 

F = Faraday (96.487 Coulomb): 


FIGURA 2 — Voltamograma cíclico para a redução 
reversível da espécie Ox (Ox e Red 
são solúveis) 
E, — potencial de inversão 
' al a 1 : ' á 

E" É ,— potencial de meio pico (catódico, anódico) 
pis Pi é 

É |,» — potencial de meia-onda 


1, 1º — corrente de pico (catódico, anódico) 
p 


(O indice o significa que as correntes são medidas em 
relação ao eixo de corrente nula) 
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R = constante universal dos gases perfeitos 
(8.315 Joule/mole grau); 

T = temperatura absoluta ("K); 

A = área do eléctrodo de trabalho (cm 9; 

C = concentração do reagente no seio da solução 
(mole/ 1); 

D = coeficiente de difusão (cm * | Ss); 

v = velocidade de varrimento de potencial (V/s). 


O potencial de pico Ep (V) para este processo, é 
independente da concentração do reagente e está re- 
lacionado com o potencial de meia-onda polarográfi- 
ca, E | 2 pela equação (1): 


Ep= E1/2— 1,11 RT/nF 


O valor de E,» pode determinar-se a partir do 
voltamograma de um processo reversível simples, 
tendo em atenção que ele ocorre num ponto corres- 
pondente a 0,85 ip (7). 

O potencial de meio pico E ,;2 é um ponto de refe- 
rência conveniente e está relacionado com E ,,2 e E 
pelas equações seguintes (7): 


p 


Ep 2= E//2+1,09 RT/nF 


AE=E,2— Ep=2,20 RT/nF 


Ox +ne  =Red 


E 


q 
a 
Ee 
“O 

Rar 
Ea! 
Q 


CORRENTE 


As equações acima só são válidas para um processo 
de eléctrodo reversível, em que no início só existe a 
espécie Ox e em que o produto Red, também é 
solúvel. 

Para as transferências de carga reversíveis, em que 
ambas as espécies Ox e Red são solúveis, de que o 
caso acabado de referir é um exemplo, Mamantov e 
colaboradores (11), referem que em voltametria de 
varrimento linear, a representação gráfica de log 
(ip— 1)/1 versus potencial, deve ser linear na gama de 
valores 0,35-0,70 ip, com um declive de 0,58 (nF/RT) (7). 

A teoria da voltametria de varrimento linear, que 
inicialmente foi desenvolvida para a análise de pro- 
cessos reversíveis simples, tem vindo a alargar o seu 
leque de aplicações, estendendo-se actualmente a 
processos de eléctrodo irreversíveis (1,4,7,12,33,34,39), 
reacções químicas acopladas a transferências de carga 
reversíveis ou irreversíveis (2,7,14-23,38,42,54), reac- 
ções químicas entre duas transferências de carga (14-16, 
26,43,45,46), transferências de multicargas (26,37,50), 
adsorção de reagentes e produtos (25,28,30,53) e di- 
versos outros estudos, dentre os quais podemos ainda 
salientar a catálise redox (19,40,42) e a determinação 
da estequiometria de complexos metálicos (47). 

Nesses estudos teóricos, critérios de diagnóstico 
apropriados, tais como a variação de ip /v !/2, Epe 
ipê / ipScom a velocidade de varrimento do potencial 
(v), têm vindo a ser estabelecidos para a caracteriza- 
ção qualitativa e quantitativa de sistemas de eléctro- 
dos desconhecidos. Uma compilação desses critérios 
é apresentada na secção seguinte, com que termina- 
mos o presente trabalho. 

A voltametria de varrimento (linear e cíclico) tornou- 
-se assim uma ferramenta importantíssima, na análise 
cinética das reacções parciais que ocorrem nos pro- 
cessos de eléctrodo complexos, no estudo de fenóme- 
nos de superficie, na caracterização de sistemas des- 
conhecidos, etc., particularmente a partir de 1973, 
data em que foi introduzida a técnica da voltametria 
integrada de varrimento (4,26,32,35,36,37). 

Os critérios incluídos neste artigo, referem-se ao 
caso dum eléctrodo plano (difusão linear). A aplica- 
ção da voltametria cíclica aos processos ocorrentes 
em eléctrodos esféricos e cilíndricos encontra-se nu- 
ma fase menos avançada, de modo que um estudo 
análogo ao apresentado, para eléctrodos com essa 
configuração, será desenvolvido mais adequadamen- 
te num futuro próximo, 

As espécies electroactivas que intervêm nos proces- 
sos de eléctrodo, indicados a seguir, consideram-se 
sempre solúveis no electrólito, excepto quando refe- 
rência em contrário, 

Além disso, para os mecanismos la XXV, pressupõe- 
se que a espécie O é a existente inicialmente em solu- 
ção, pelo que, alguns dos critérios apresentados, res- 
peitam especificamente à redução dessa espécie O. 


CRITÉRIOS DE DIAGNÓSTICO 
Mecanismo I 


Transferência de carga reversível (1,4,7) 
O+ ne =R 
a) ip= 0,452 (nF)3/2 (Dv) 1/2 AC / (RT)!/2 


b) Ep = f(v) 
Ep Hf(C) 
c) Ep — Ep =-—(56/n) (T/298) mV 
d) ip / vi/ZAF(v) 
e) db [ ip =| 
) ES -Ep2=-2,20 RT/nF mV 


Mecanismo IH 


Tranferência de carga reversível (1,11,41,50) 
Espécie R insolúvel 


O+ne =R 


a) ip = 0,610 (nF)3/2(Dv) !/2 AC/(RT)!/2 
b) Ep * f(v) 
e) ip/ vi/2f(v) 
d) Ep=f(C) 

Ep será desviado no sentido catódico, com o de- 
créscimo da concentração da espécie O. 
O is Lip > 
O Ef —Ep/2=—0,77 RT/nE myV 
g) E é uma função linear de log (ip —1) na gama de 
correntes de 0,5 ipa 0,9 1p. 

O coeficiente angular da recta correspondente é 

dado por: 0,454 nF/RT 


Mecanismo HI 
Transferência de carga múltipla (21) 


mO +ne = gqR 


i) Reacção de 2.º ordem m=2 q-1 


a) ip = 0,353 (nF)3/2 (Dv) 1/2 AC/(RT)!/2 
b)ip / ip = 1,09 

c)ip/ vi/2 f(v) 

d) Ep f(v) 

e) Es — Epp =—(81,4/n) (T/298) mV 
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ii) Reacção de 3.ºordem m=3 q-1 


a) ip = 0,303 (nF)3/2 (Dv) 1/2 AC/ (RT) 112 
bip /ip =1,16 

Cc) ip/ vl/2 5 f(v) 

d) Ep »* f(v) 

e) EF — Ep;2= —(104,5/n) (T/298) mV 
Mecanismo IV 

Transferência de carga quase reversível (12,13) 


lento 
Õ +ne = R 
a) E p = f(v) 
b) Ep — E; =-—(60/n) (T/298) mV para baixos v, mas 
aumenta à medida que v aumenta. 
c) ip / vl/2=f(v) 
d)ip /ip =lapenaspara a=0,5 


a — coeficiente de transferência de carga. 


Mecanismo V 


Transferência de carga irreversível (1,7,12,33,34) 


o 
O + ne = im 
a)ip=0,227nFACK,exp[(— an, F/RT)(Ep -E )] 
K =Kçexp[(—- ana F/RT) (E - E?) ] 
Eº — potencial formal de eléctrodo 
K — constante de velocidade da reacção 
K« — constante de velocidade da reacção ao potencial 
de inversão 
na —n. de electrões postos em jogo no processo parcial 


que controla a velocidade 


b) E, e Ep/2 diminuem de (30 / a na) (T/298) mV 
cada vez que v decuplica 
C)ip/ vl!2xf(v) 
d)ip =0 
e) Ep — Ep/2= 1,857 (RT / ana F) 
9) Ep —Eº=-—(RT / anaF) 
[0,78 + In (Db) !'2 — InkKç] 
b=anaFv/RT 
g) (Ep)2 — (Ep) =(RT/ ana F)In(vi/ v9)!!? 
(Ep) 1 — potencial de pico à velocidade v1, 
(E p)2 — potencial de pico à velocidade va, 
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Mecanismo VI 
Transferência de carga irreversível (1,12,33) 


Espécie R insolúvel 


. K, 
Ô + ne — R 


a)ip=0,5nF%2 (a n)!/2 (Dv)!/2 AC/(RT)!!2 
b) E, e Ep/2 diminuem de (30 / a na) (T/298) mV 


cada vez que v decuplica. 


Mecanismo VIH 


Transferência de carga reversível seguida 
duma reacção quimica reversível (7,4) 
O+ne =R 


Kr 
R = Z 
Kb 
kK =Kr/ Kb Kp+Kç=1 a=nFv/RT 


Considerando o caso (*) em que le K são elevados 


a) ip = 0,878 (nF)3/2 (Dv) 1/2 AC/(RT)!/2 

b) Ep desloca-se no sentido catódico à medida que v 
aumenta, e esse desvio é de cerca de (60/n) (T/298) mYV, 
cada vez que v centuplica. 

c)ip/ vl/2 5 f(v) virtualmente ** 

d) a / is decresce a partir de | à medida que v aumenta 
e) Ep, — E,/2 aumenta quando v aumenta 


Mecanismo VHI 


Transferência de carga reversível, seguida de 
dimerização reversível (20) 
O+ne =R 


* Para além deste, é possível considerar neste meca- 
nismo os dois casos seguintes, de acordo com o valor 
do parâmetro cinético Kal/2 / 81/2 

1) Se a reacção quimica for muito rápida, a onda 
deslocar-se-d no sentido anódico, ao longo do eixo 
dos potenciais. 

2) Sea reacção quimica for muito lenta, (Kj >> kb). 
pode considerar-se irreversível (Mecanismo X). 

** Esta expressão significa que o critério não é aplicável a toda a 
gama de valores possíveis de v ou c. 


i) Se K é pequeno (corrente cinética pura) 

a) ip = 0,527 (nF)3/2 (Dv) I/2 AC/(RT)!/2 

b)ip/ vl/2f(v) 

c) Ep = f(v), variando de (20/n)(T/ 298) mV no senti- 
do catódico cada vez que v decuplica 

d) Ep =f(C) variando de (20/n) (T/298) mV no senti- 


do anódico cada vez que € decuplica 


ii) Se K é elevado (corrente de difusão pura) 

a) ip = 0,50 (nF)3/2(Dv) 1/2 AC/(RT) 1/2 

b) Ep =f(v) 

c) Ep =f(C) variando de (30/n) (T/298) mV no senti- 
do anódico, cada vez que €C decuplica 


Mecanismo IX 


Transferência de carga reversível, seguida 
de reacção multimolecular reversível (7,17) 


O+ne =R 


Kr 
mR = qZ 

Kb 
m=3 m=2 m=1 m=2 m=1 m=1 
q=1 | q=l o q=2 q=3 


a) Para ip / io independente de v, tem-se, consoante 
os valores de m e q, os seguintes valores de é / a 


respectivamente, pela ordem indicada: 

dg | ip = 1,143; 1,087; 1,0; 0,952; 0,917; 0,869 
Geralmente, verifica-se um aumento de ip, à medida 
que q/m diminui. 

b) Os valores de (Ep — E,,2) em mV, serão pela or- 
dem indicada, os seguintes: 
(—35/n)(T/298); (—40,2/n)(T/298); (-56/n)(T/298); 
(—70,3/n) (1/298); (—83,9/n) (T/298); 

(—108,6/n) (T/298) 

Em geral verifica-se um decréscimo de Ep — Ep/2 à 


medida que q/m diminui, 


Mecanismo X 


Transferência de carga reversível, seguida duma 
reacção química irreversível (7,18,22,23) 


O+ne =R 


REZ 
Se Kr/ a é pequeno, o voltamograma que se obtém 
é o característico duma reacção de eléctrodo reversi- 
vel, uma vez que é desprezável o efeito da reacção 


química. 
a=nFv/RT 


Se Kr/ a é grande, não se observa corrente anódi- 
ca, tal como acontece no caso da transferência de 
carga totalmente irreversível. 


a) Para baixos valores de v, E, torna-se mais catódi- 


co, diminuindo de cerca de (30/n) (T/298) mV cada 


vez que v decuplica. 


A referida variação é menos acentuada à medida que 
v aumenta. 

b) Ep = f(C) 

c) ip / v!/2 2 f(v) virtualmente 


d)ip / ip tende para | quando v aumenta. 


Mecanismo XI 


Transferência de carga reversível seguida duma 
reacção de dimerização irreversível (7,18,22,23) 


O+ne” = R 


Kr 
2 2 


C " ã Ai 
a) Ep torna-se mais catódico à medida que v aumenta, 


podendo essa variação traduzir-se por: 
AEp /Alogv=(20/n) (T/298)mV/unidade de log v 
b) E torna-se mais anódico à medida que €C aumen- 
ta, podendo essa variação traduzir-se por: 
PA ES [Alog C =(10/n)(T/298) mV/unidade de log € 
)ip/ v!/2xf(v) 
d) EF — Ep/2= —(39/n) (T/298) mV 
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Mecanismo XII 


Transferência de carga reversível seguida de 
reacção cataliítica irreversível (7) 
O+ne” = R 
Ko 
R+Z — OQ 
a) Ep torna-se mais catódico, sendo essa variação de 
(60/n) (T/298) mV cada vez que K, / a decuplica 
b) Ep2*f(v) seKç/a > 10 
c) O pico alarga-se à medida que Kç / a aumenta, 
sendo impossivel detectá-lo para Kç /a> 1. 
d) ip / v!'2 decresce à medida que v aumenta, para 
baixos valores de v e é independente de v para eleva- 


dos valores de v. 
e) / sê = | desde que E, < E. (35/n) mV 


onde E, — potencial de inversão 


Mecanismo XII 


Transferência de carga irreversível seguida duma 
reacção catalítica irreversível (7) 


e 
O+ne — R 


a) Ep e Ep/2 tornaram-se mais catódicos à medida 
que Kc/ b aumenta. 


b) ip / v!/2 decresce à medida que v aumenta, para 
baixos valores de v 


ip / v!/2 &* f(v) para elevados valores de v 


c) Não se observa corrente anódica. 


Mecanismo XIV 


Reacção ECE com dismutação do produto da pri- 
meira transferência de carga (15) 


O+nie ==R 
Kr 
R == A 
Kp 
Ka 
IR — 0+4Z 
A +n9€ = 


a) À variação de ip com a concentração inicial C 
permite distinguir o mecanismo ECE simples, daque- 
le em que há dismutação. 


No caso ECE simples, 1, a C€, enquanto que quan- 
do há dismutação acoplada, as variações de ipcom €, 
são maiores que o que seria de esperar duma propor- 
cionalidade directa. 


b) A variação de ip com v neste caso, é diferente da 
que se verifica no mecanismo ECE simples. 

Fazendo variar ve Ce representando ip em função 
do parâmetro cinético À = f(1/v), é possível detectar a 
etapa determinante da velocidade e determinar a sua 
constante de velocidade (K 4), no caso da dismutação, 
ouK'/ K” no caso do mecanismo ECE simples (onde K' 
representa a constante de velocidade do processo 
globale K =Kr/ Kbp). 

Finalmente, no caso do mecanismo ECE com dis- 
mutação, verifica-se que: 


Kg * f(v) virtualmente 
Kg f(C) virtualmente 


Mecanismo XV 


Transferência de carga reversível com 
adsorção do reagente e do produto (25) 


O = Oads 
O+ne” = R 
K. R ads 


Se o reagente ou produto é fracamente adsorvido, 
verifica-se um acréscimo das correntes de pico, o ê 
Ds em comparação com as obtidas para a transfe- 
rência de carga reversível simples sem adsorção (Me- 
canismo 1). 

Se o reagente ou produto é fortemente adsorvido, 
poderão ocorrer picos de adsorção (iadgs) bem dife- 
renciados depois ou antes do pico de redução ou de 
oxidação. 


a-l) Reagente fracamente adsorvido 


Neste caso, ip aumenta quando v aumenta e dimi- 
nui quando C aumenta; ip / ip é aproximadamente 
constante e geralmente diferente de 1. 

Para velocidades de varrimento baixas, contudo, a 
influência da adsorção na corrente torna-se desprezi- 
vel e obtém-se os mesmos resultados que para o caso 
reversível em que is [ E = 1, 

Quanto à influência da concentração, verifica-se 
que a fracção da corrente de pico, devida à substância 
adsorvida, decresce com o aumento de concentração. 


a-2) Produto fracamente adsorvido 


Nestas condições, io aumenta quando v aumenta, 
conservando-se ip praticamente constante. 


b-1) Reagente fortemente adsorvido 


Este caso é semelhante ao anterior, com a diferença 
de que surge um pico de adsorção que, no ramo cató- 
dico, aparece após o pico de difusão, tal como se 
ilustra na Fig. 3-C, na página seguinte. 

Por este facto, abstemo-nos de quaisquer conside- 
rações sobre o assunto. 


b-2) Produto fortemente adsorvido 


Neste caso, o comportamento do sistema para baixos 
valores de v aproxima-se duma redução reversível 
normal. À medida que v aumenta, aparece e vai 
aumentando iags, € ip Vai decrescendo, de tal modo 
que, para valores de v elevados, 1p tem um valor des- 
prezável (correspondendo ao desaparecimento do pi- 
co de difusão), e só o pico de adsorção lag; acaba por 
ser visível. 

Quando C aumenta, ads diminui ocorrendo a um 
potencial mais anódico, enquanto que ip só varia 
ligeiramente. 


Os voltamogramas cuja simetria aumenta quando 
v aumenta ou € diminui e que, em geral, apresentam 
um 2.º pico, localizado imediatamente antes (pré- 
-pico) ou imediatamente a seguir (pós-pico) ao pico 
de difusão, são indício de adsorção. 

Se os processos de adsorção e difusão ocorrem si- 
multaneamente não se observam tais picos. 

Na presença de adsorção é também costume observar- 
-se um aumento de 1p / € quando € diminui. 

Contudo, tanto para elevados €, como para baixos 
v, a influência da adsorção nos voltamogramas é rela- 
tivamente desprezível. 


FIGURA 3 — Voltamogramas teóricos para o mecanismo 
de transferência de carga com adsorção (25) 


— Reagente fracamente adsorvido 
— Produto fracamente adsorvido 


Reagente fortemente adsorvido 


vor > 
| 


— Produto fortemente adsorvido 


(As linhas a tracejado referem-se a uma transferência de 
carga reversivel sem adsorção) 


Mecanismo XVI 


Transferência reversível de multicargas 
com regeneração catalítica (R-C-R) (42) 


A+ne” = B+Z Ef 
Kc 
B+ne = €C E» 
nj=n92=n 


a) Nos casos em que a onda correspondente à 1. 
transferência de carga, é bem separada da correspon- 
dente à 2.º transferência de carga, 
(/AEº/>(180/n)(T/2988)mVcoma Eº=E3 - Ep <0) 
verifica-se que aquela apresenta as características de 
uma reacção catalítica simples, na ausência da 2: 
transferência de carga. 


a-l) Neste caso, se Kr / a é pequeno, a reacção 
quimica não tem qualquer efeito no mecanismo, 
obtendo-se um voltamograma idêntico ao de uma 
transferência de carga reversível. 

À medida que Ky / a aumenta, o efeito da reacção 
química torna-se cada vez mais importante e a cor- 
rente aumenta (devido à regeneração química do 
reagente). 


e] 
a 
[ee 
om | 
[ui 
tu 
Õ 
Q 
q | 
q 
Q 
103 
t+ 
[e 
= 
EL 


0,1 0,0 -O,1 


Potencial, mV 
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Para valores de Kr / a elevados, não se observa 
qualquer pico. O polarograma apresenta um patamar 
cuja corrente é proporcional a Ko e independente 
de v. 


a-2) Ep (ou Ep/2 quando não se observa pico), 
desvia-se no sentido catódico de cerca de (60/n) 
(T/298) mV cada vez que K, / a decuplica. 


a-3) Se o potencial de inversão, (E,), se encontra 
entre os das duas ondas, e é pelo menos (35/n) 
(T/298) mV mais catódico que o potencial de pico da 
onda 1, (Epi), o voltamograma obtido, é idêntico ao 
de uma reacção catalítica acoplada a uma única 
transferência de carga reversivel. 

— Se E, é (65/n) (T/298) mV mais catódico que Epa, 
aparece um pico anódico (ip2) correspondente à 2? 
transferência de carga, cuja altura decresce à medida 
que K« / a aumenta. 


b) Se a sepração entre as duas ondas é pequena 
(/ A Eº/<(180/n)(T/298) mV com À Eº< 0): 

— Para pequenos valores de Kr / a, as duas ondas 
misturam-se distorcendo-se. 

— À medida que Ky / a aumenta, desaparece qual- 
quer vestígio de diferenciação entre as duas ondas, 
aparecendo em sua substituição um patamar. Este 
patamar é contudo mais estreito que o referido em a-1). 


c) SeA Eº=00u A Eº>0Q 
Não se observa efeito catalítico. 


Aparece somente uma onda catódica e outra anódica. 


d) A forma e a posição das ondas variam com nj € no. 


Mecanismo XVII 


Transferência de multicargas com regeneração 
catalítica, em que apenas a |“ transferência 
de carga é reversível (R-C-1) (42) 


Arne” = B+Z 


Kc 
B + me = E 
nj=n92=n 
— À separação entre a onda correspondente à 1. 


transferência de carga e a referente à 2.º transferência 
de carga depende de K, / a. 


a) Se o pico correspondente à 1.º transferência de 
carga é bem distinto do correspondente à 2º transfe- 
rência, isto é, se A Eº<Qe 
[A Eº/>(180/n)(T/298) mV 


ER à 
a-l) 1p2 é menor que no caso do mecanismo XVI. 
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a-2) Para valores baixos de K / a(Kç/ a< 0,1), 
observam-se 2 picos catódicos distintos. 


a-3) Para valores mais elevados de Kç / a 
(K./a> 0,1), ip2 é muito baixo enquanto a onda | 
apresenta um largo patamar, semelhante ao observa- 
do para o mecanismo XVI. 

a-4) Se / Epp />/ Es/>/ Epil OS voltamogra- 
mas obtidos são idênticos aos do Mec. XVI (a-3). 

Se / Eç/>/ Epa / — não se observa resposta anódica. 


b)SeA Eº<0e/AEº/<(180/n)(T/298) mV 


b-1) Para valores baixos de Kc« / a(Kc/a<0,1),as 
ondas sobrepõem-se e aparecem distorcidas. 

b-2) Para valores mais elevados de Kc /a(Kç/ a 
> 0,1),a sobreposição das ondas origina um patamar 
e E p2 desvia-se de (30 / a na) (T/298) mV no sentido 
catódico, cada vez que v decuplica. 

c) Se a separação entre as duas ondas é nula 
A Eº=0,ou positiva 4 Eº > 0, os voltamogramas 
obtidos são idênticos aos do Mec. XX. 


Mecanismo XVII 


Transferências de carga reversíveis 
consecutivas (R-R) (50) 


One = Ri Ey 
Ri+ne = R> E; 


aJSeA Eº<Q0e/A Eº/>(118/n)(T/298) mv, 
onde n = ny = no (Fig. 4-4) e se O for uma espécie 
mais facilmente redutível que R|, o voltamograma 
cíclico apresenta 4 ondas separadas, (2 corresponden- 
tes à redução global de 0a R$€e 2 à oxidação de R$ a 0). 

a-1) Ip [ viZ & f(v) Critério característico de 
ausência de reacções químicas acopladas às reacções 
de transferência de carga. 


a-2) En2 -— Epi * f(v) 


fig É E : E otã A 
a-3) ipi / ipi = ip2 / ip? = ol [àp2 =ipl / àp2 =] 


a-4) Desde que a altura do pico de cada onda seja 
medida em relação a uma linha de base determinada 
pelo ramo descendente da onda que a precede: 


ip2 ) ipi A t (Ep> Na Epi) 
a-5) ipa (n))'2 


b)SeAEº<0Qe/AEº/< (100/n)(T/298)mV 
(Fig. 4-B) as ondas correspondentes às duas transfe- 
rências de carga aparecem fundidas numa única onda 
distorcida, cuja altura e forma não são características 
de uma reacção reversível. Esta onda distorcida asseme- 
lha-se a uma onda irreversível, da qual contudo esta 


FIGURA 4 — Voltamogramas cíclicos para diferentes 
casos do mecanismo R-R (50) 


no/nj=1 

A — A. Eº = — 180mV 
B— A Eº= -9mV 
C— A Eº = O mV 


D— A Eº=+ 180mV 


onda se distingue pelo facto de não variar a sua posi- 
ção no eixo dos potenciais em função de v. 


c) Se A Eº=0(Fig. 4-C), a onda respectiva apre- 
senta um só pico anódico e outro catódico de altura 
intermédia entre a referente a uma transferência de 
carga de | electrão e a correspondente a 2 electrões 
(no caso den, =n2 = 1) (ip para uma transferência de 
2 electrões é (2)*'/2 vezes maior que lp para uma trans- 
ferência de 1 electrão, no caso de ambas as transfe- 
rências serem reversíveis). A onda resultante da so- 
breposição das duas transferências de carga é caracteri- 
zada por um Eº=(E, +E; )/2. 

d) Se À Eº> 0 aparece um só pico anódico e | outro 
catódico. (Fig. 4-D). 


Mecanismo XIX 


Transferências de carga irreversíveis 
consecutivas (1-1) (50) 


O +nj; E = R, 
Rj+n9 e” doi Ro 


a) Se as duas espécies 0 e R, forem reduzidas a 
potenciais suficientemente diferentes, cada uma das 
transferências de carga apresenta as características 
duma transferência de carga irreversível isolada (Me- 
canismo V). 

Não há resposta anódica, apresentando o voltamo- 
grama 2 picos catódicos distintos. 


b) Se o potencial de redução da espécie 0, coincidir 
com o de R,,0 voltamograma apresenta uma única on- 
da, com as características resultantes da combinação 
das duas transferências de carga envolvidas. 

Assim, surge apenas um pico catódico, verificando- 
se a ausência de picos anódicos. 


Tanto no caso a) como em b) verifica-se que: 


ip / v!/25 f(v) qualquer que seja o valor de v 


Mecanismo XX 


Transferências de carga consecutivas 
em que só a |./ é reversível (R-1) (50) 


O+ne” = Rj 


Rj+n9e = R> 


a 
vim 
[ 
e] 
[a 
te 
Q 
Q 
q 
a 
Q 
te] 
E» 
E 
pes | 
E 


-200 200 


Potencial, mV 


a) Se os potenciais de redução das espécies 0 e R;, 
forem suficientemente diferentes, observam-se duas 
ondas no voltamograma, que podem ser consideradas 
como resultantes das duas reacções separadas. 


a-l) A 1 transferência de carga tem as caracteris- 
ticas já referidas para o Mecanismo 1, enquanto a 2; 
transferência apresenta as características do Meca- 
nismo V. 

Deste modo, o voltamograma obtido apresenta 2 
picos catódicos e um só pico anódico, corresponden- 
do este à 1.º transferência de carga. 

a-2) ip / v!/2 f(v) 

b) Se os potenciais de redução de 0 e R; forem 
coincidentes observa-se uma única onda resultante da 
combinação das duas ondas (correspondentes a cada 
uma das transferências de carga). 


b-1) O voltamograma apresenta um só pico anódi- 
co e um só pico catódico. 
b-2) ip / vl!2 f(v) 


Mecanismo XXI 


Transferências de carga consecutivas, 
em que só a 2.º é reversível (I-R) (50) 
O+ne' —- R; 

Ri+ne = R5 


a) Se os potenciais de redução das espécies 0 e R; 
forem suficientemente diferentes, o voltamograma 
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correspondente a este mecanismo apresenta duas on- 


das que podem ser consideradas como resultantes das 
duas reacções de transferência de carga separadas. 


a-1) A 1.º transferência de carga apresenta as carac- 
terísticas referidas para o mecanismo V, enquanto a 
2º transferência se identifica com o mecanismo 1. 

O voltamograma obtido apresenta 2 picos catódi- 
cos e um só pico anódico, este último correspondente 
à 2º transferência de carga. 


a-2) ip / vl/2 f(v) 
b) Se os potenciais de redução de 0 e R;, forem 


coincidentes observa-se uma única onda, resultante 
da combinação das duas ondas referidas. 


b-l) O voltamograma apresenta um único pico 
anódico e um único pico catódico. 


b-2) ip /V 2% f(v) 


Mecanismo XXII 


Reacção química irreversível entre duas 
transferências de carga reversíveis 
(7,14) 


OU+nje = KR E, 


R = C 


C+ne” = D E; 


| ns , E E 1 3/2 
a) Se n2 / nj = constante, ipja ny eipza nz”. 
E | 
lpl — corrente de pico catódico para a 1.º transferência 
de carga (0 — R) 


.€ atua 
Ip2 — corrente de pico catódico para a 2.º transferência 
de carga (C — D) 


b) O valor relativo dos dois picos catódicos pode 
utilizar-se na determinação de Kg recorrendo-se a 
uma curva 


Jc E 
Ipi/ 1p2 em função de Kr-/ a 


c)ip / v!/2 = f(v). Este critério pode usar-se para 
distinguir este mecanismo da transferência de multi- 


cargas. 
d) Se A Eº<O0 


d-1) Para pequenos valores de Kç / a, o primeiro 
processo de transferência de carga tem as proprieda- 
des da transferência de carga reversível, seguida du- 
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ma reacção química irreversível. Assim, só aparecem 
bem evidenciados um único pico catódico e um só 
pico anódico, correspondentes ao 1.º processo de 
transferência de carga. 


É , €£ E ii 
Lpi/ vl2 =f(v) virtualmente; | p3/ vi/2 =f(v) =0 


ip2 [ vtZ =f(v) =0 ipi | vll2=f(v) 


d-2) Para elevados valores de Kç / a, observam-se 
dois picos catódicos bem acentuados; a resposta anó- 
dica agora é simplesmente a correspondente à segun- 
da transferência de carga reversível. 


a 
ip2/ vl/2=f(v) HO 


.d 


ip] [w 1/2 =f(v); ip? Ev 1/2 =f(v) 


e) Se A Eº=0 


— Qualquer que seja o valor de Kç / a, verifica-se a 
existência dum único pico catódico e dum único pico 
anódico. A altura do pico catódico diminui com o 
aumento de v, enquanto a do pico anódico se mantém 
praticamente independente de v. 

N)SeA Eº>O 

f-1) Para pequenos valores de Kç / a, aparecem 2 
picos catódicos e somente um pico anódico corres- 
pondente à 1.º transferência de carga, uma vez que o 
pico anódico corresponde à 2.º transferência de carga, 
vai diminuindo à medida que Kç / a diminui desapa- 
recendo praticamente para valores muito baixos de 
Kç/ a. 


ipi / vH/2=f(v) ip2/ vl/2=f(v) 


Ao, mw 
ipi / v!/2=f(v)£0 Lpz / vZ =t(v) = 0 


f-2) Para elevados valores de Ky / a aparecem 2 
picos catódicos. Contudo, só aparece um único pico 
anódico correspondente à 2.“ transferência de carga. 


ipi / vl/2 = f(v) ipz / v!/2 = f(v) 
ipi =0 ip2 / v//25 f(v) 
Mecanismo XXIII 


Reacção quimica irreversível, precedida 
por uma transferência de carga reversível 


e seguida duma transferência de carga 


irreversível (7,14) 


U+ne” = R E, 


A forma dos voltamogramas obtidos para este me- 
canismo depende da separação entre duas ondas, do 
número de electrões envolvidos, e do parâmetro ciné- 
tico Kç/ a. 

Além disso, as características da 2.º transferência 
de carga, dependem de a n;/ nj. 


a)SeA Eº<OQ 


a-l) Para valores de K; / a muito baixos, o volta- 
mograma apresenta praticamente só um pico anódico 
e um catódico, correspondentes à 1“ transferência de 
carga, visto que ambos os picos (catódico e anódico) 
correspondentes à 2. transferência de carga, vão di- 
minuindo à medida que v aumenta, até que a partir 
de determinado valor de v deixam de aparecer: 


ipi [via f(v) (virtualmente) ip2 / vi? =f(v)=0 


ipi /v/2=((0H0 ibo/vi2=Hv=0 


a-2) Para valores de Kr / a muito elevados apare- 
cem somente 2 picos catódicos, ipi € ipa, 

Ausência de picos anódicos. 
€ | : e 
ipi / vI/2 f(v) (virtualmente) — ipz/ vi/2 =f(v) 


a-3) Ep2 desvia-se no sentido catódico à medida 
que v aumenta, dependendo este desvio do valor de 
Kç/a. 


b)SeA Eº=0 
Neste caso as duas ondas sobrepõôem-se, aparecen- 
do assim um único pico anódico e um único pico 
catódico para valores baixos de Kç / a. 
oc 
ip / vl/2= f(v) 
a 
Lp / vl/2 = f(v) 


Para valores elevados de Kr / a obtém-se somente | 
pico catódico: 


a 


c)Se A Eº>0 


Este caso é análogo àquele em que À Eº< 0 com as 
seguintes alterações: 


E E | ; 
A variação de ipi / v!/2 com v é agora mais acentuada. 


* ma Re Í 
A variação de ip> / v!/2 com v é agora menos acentuada. 


Mecanismo XXIV 


Reacção quimica irreversível, precedida 
por uma transferência de carga irreversível 
e seguida duma transferência de 

carga reversivel (7,14) 


0O+ne — R 
Kr 
R = C 


C+ne = D 


A 1º transferência de carga não é afectada pela 
existência da reacção química e apresenta caracteris- 
ticas idênticas às do mecanismo YV, 

A 2º transferência de carga é semelhante à homó- 
loga do mecanismo XXII. 


a)Se A Eº<Q 

— Em toda a gama dos possiveis valores de v, 
ipi f(v) eipi=0. 

a-l) Para valores baixos de Ki / a, aparece só um 
pico catódico correspondente à 1.º transferência de 
carga, uma vez que tanto o pico catódico como o 


anódico, correspondentes à 2.º transferência de carga, 
vão diminuindo à medida que v aumenta. 


ipa / vl/2=f(v)= 0 
ip2 / vl/2=f(v)= 0 


a-2) Para valores elevados de Kg / a, surgem dois 
picos catódicos (um correspondente à 1.º e outro à 2º 
transferência de carga), e um só pico anódico, corres- 
pondente à 2.º transferência de carga. 


ip2 [vliZ=f(v) e i p2 [vl2=f(v) 


b) Se A Eº=0 


— Neste caso há sobreposição das ondas corres- 
pondentes a cada processo de transferência de carga. 


b-1) Para valores baixos de kr / a, só aparece um 
pico catódico, uma vez que o pico anódico vai dimi- 
nuindo à medida que v aumenta: 


ip= (WO ip / vi/2=f(vy)= 0 


b-2) Para valores elevados de kr / a, aparece um 
pico anódico e outro catódico: 
K + a 
Ip=f(vV) HO e Ip=[(V)ZO 
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c)SeA Eº>O0 


Este caso é análogo ao caso a) (A Eº< 0), comas 
seguintes diferenças: 


ipi [v 12 = f(v) 


e : 
Ip2 assume valores menores que a homóloga no 
caso a), em toda a gama de valores de v. 


Mecanismo XXV 


Reacção quimica irreversível, 
precedida e seguida duma transferência 
de carga irreversível (7,14) 


O+nj;e” — " R 


Cr+ne” > D 


— O voltamograma é análogo ao resultante da 
combinação da onda correspondente à 1.“ transferên- 
cia de carga do mecanismo XXIV, com a correspon- 
dente à 2.º transferência de carga do mecanismo XXIII. 

A 1º transferência de carga não é influenciada pela 
reacção química. 

A 2º transferência de carga depende de kr / a 


a)Se A Eº<0 

Não há resposta anódica, qualquer que seja o valor 
de v. Para elevados valores de kr / a aparecem 2 picos 
catódicos. Contudo, enquanto que 
e e o doe a 
ipi / vi/25 f(v); ip2 / v!/2 diminui à medida que v 
aumenta, assumindo valores muito baicos para ele- 
vados valores de v. 


ip1=0 ip2=0 
b) Se A Eº=0 


Observa-se um único pico catódico, não havendo 
qualquer resposta anódica, qualquer que seja o valor 
de v. 


ip/vi2=HWAOo ip=0 


c)Se A Eº>0 


Não ha resposta anódica, qualquer que seja o 
valor de v. 

Para elevados valores de kr / a, só se observa um 
pico catódico, correspondente à 1.“ transferência de 


carga (ipi ). 
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. p E 
Para baixos valores de kr / a, surge além de i pj, um 
er pao ? 
outro pico catódico (1 p2), correspondente à 2. trans- 


ferência de carga. 


ipi / v!2=f(v) ipa) vi/2=f(y) 


a a 
Ip1=0 1p2=0 


Mecanismo XXVI 


Reacção quimica de |“ ordem reversível, seguida 
duma transferência de carga reversível 


(4,7,20,38,48) 
Kr 
A = O 
Kb 


a) Ep =f(v) para elevados valores de v, desviando-se 
de (60/n) (T/298) mV no sentido anódico, cada vez 
que v centuplica. 


E, é praticamente independente de v, para baixos 
valores de v. 


b)ip/v !/2 decresce quando v aumenta. 


| JE 
c) ip / 1p> | e aumenta quando v aumenta. 
Ip/ip=1 para baixos valores de v. 


d) ip = f(v) ip decresce à medida que v aumenta. 


Mecanismo XXVII 


Reacção de monomerização reversível, 
seguida de transferência de carga 
reversível (20) 


i) Se K é elevado (corrente de difusão pura) 
a) ip = 1,087 (nF)3/2 (D,v)!/2 AC, (RT)? 
b) Ep 2 f(v) e Ep* f(v) 
)E2=N€) é Ep=HC;) 
d) ip / vl/2 f(v) 


o Os eo O ———O 


ii) Se K é pequeno (corrente cinética pura) 


a) ip = 1,155 nAF (D,Kp) K34 C)!/ (RT)!!? 

b) Ep;2 varia de (30/n) (T/298) mV no sentido 
anódico, cada vez que v decuplica. 

c) Ep/2 varia de (15/n) (T/298) mV no sentido ca- 


tódico, cada vez que €, decuplica. 


Mecanismo XXVIII 


Reacção química reversível, seguida 
duma transferência de carga 
irreversível (7) 


Kr 
Z - 0 
Kb 
K' 
O+ne' — R 


a) Não existe resposta anódica. 


b) ip / v!/2 decresce quando v aumenta. 

c) Para valores elevados do parâmetro cinético 
b!'2 | K11'2, Ep desloca-se no sentido anódico de 
cerca de (60 / a na) (T/298) mV cada vez que v 
centuplica. 

E, é independente do parâmetro cinético, quando 
este assume valores baixos. 


d) Ep/2 aumenta quando v aumenta. 
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SUMÁRIO 


O artigo refere a estrutura organizacional do Curso 
de Engenharia de Minas do IST, mencionando os 
seus principais laboratórios e centros de investigação. 
São apresentadas as principais linhas de I&D em rea- 
lização nesses laboratórios e centros de investigação, 
bem como os principais casos de colaboração com o 
exterior (laboratórios nacionais, administração pu- 
blica e indústria). 


ABSTRACT 


The organizational structure and R&D activity of 
Mining Department of IST are presented. Some 
comments about the mining engineering course are 
also made. The main R&D activities of the Mining 
Department laboratories (Laboratory of Mineralogy, 
Laboratory of Geology and Laboratory of Ore Dres- 
sing are emphasized. A short presentation on the 
main research projects undertaken in those laborato- 
ries is made. 


A endogenização criteriosa e crescente das bases 
científicas e tecnológicas que suportam a Civilização 
actual é condição necessária para o nosso desenvol- 
vimento. Entenda-se por desenvolvimento a evolução 
global de uma colectividade para uma direcção posi- 
tiva o que pressupõe formas de modificação política, 
econômica, social e técnica com a satisfação o mais 
equitativa das necessidades físicas, psíquicas, cultu- 
rais e morais no geral e em particular no referente à 
autonomia de decisão e de gestão, à alimentação 
equilibrada, ao emprego seguro e convenientemente 
retribuído, à habitação conforme com as normas cul- 
turais, a nível de saúde correcto e à liberdade de 
expressão. 

É imperioso que se acelere o processo de endogeni- 
zação técnico-científico entre nós. O grande proble- 
ma com que nos defrontamos traduz-se, fundamen- 
talmente, na falta de relação entre as actividades 
devotadas à investigação científica, ou seja à pera- 
ção de conhecimento e a evolução das técnicas de 
produção que acontece arrastarem-se rotineiramente 
na indústria do país. 

Eis aqui um vasto campo a colmatar. Uma Univer- 
sidade Técnica e muito particularmente o Instituto 


Superior Técnico pelos seus docentes e investigadores 
deve procurar, com entusiasmo e sem perda de tem- 
po, fazer a ponte que permita o aumento da inovação 
tecnológica, 

É bom ter presente que o conhecimento científico e 
técnico disponível em um pais e o modo como é usa- 
da nos processos de produção são factores de poder 
com forte incidência nas relações desse país com os 
restantes. Eis porque, nas relações políticas entre os 
estados, são factores fundamentais de exercício do 
poder o grau de superioridade ou de dependência 
cientifico-tecnológica com que eles dialogam. 

Do exposto se conclui que em uma escola de enge- 
nharia, a ID&D aí em curso deverá procurar dar 
resposta aos problemas em equação na sociedade em 
que se insere. Para isso deverão ser seleccionados 
problemas a resolver tendo em consideração a possi- 
bilidade de obter resultados concretos em um periodo 
de tempo razoável expresso em termos de produção 
de tecnologia utilizável em tempo útil pela comunidade. 

Na sequência das ideias expostas, os laboratórios 
que suportam o Curso de Engenharia de Minas do 
IST conduzem alguma investigação cujas principais 
linhas vamos apontar a traço grosso. 


A organização do Curso de Minas do IST, a traços 
largos é a seguinte: há três laboratórios: o Laborató- 
rio de Mineralogia e Petrologia (LAMPIST), o La- 
boratório de Geologia (LAGIST) e o Laboratório de 
Mineralurgia e Planeamento Mineiro (LMPMIST) 
que suportam a maioria das cadeiras do curso dando- 
-lhe a base laboratorial e prática. O Grupo de Explo- 
ração de Minas agrega as disciplinas de Lavra de 
Minas e Geomecânica. 

O Laboratório de Mineralogia e Petrologia do IST- 
-LAMPIST — está razoavelmente equipado para 
trabalhos de microscopia em luz transmitida e em luz 
incidente e para análises físicas e químicas de mine- 
rais, minérios e rochas. Com efeito possui aparelha- 
gem de análise térmica diferencial e ponderal, de di- 
fracção pelos raios-X, de fluorescência pelos raios-X, 
de espectrografia pelos raios infra-vermelhos e de es- 
petrografia de absorção atómica. Tem algumas faci- 
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lidades de cálculo que lhe são facultados por um mini- 
-computador próprio. No LAMPIST está sediado o 
Museu Bensaúde que possui importantes colecções de 
minerais e rochas nacionais, além de outras estrangeiras. 

Junto do LAMPIST funciona o Centro de Petro- 
logia e Geoquímica da UTL suportado financeira- 
mente pelo INIC. 


O Laboratório de Geologia (LAGIST) do IST 
dispõe de equipamento microscópio para estudo de 
minérios e de algum equipamento para prospecção 
geofisica (eléctrica, sísmica e gravimétrica) e de uma 
instalação para estudo experimental de métodos eléc- 
tricos. Possui um museu com colecções de estratigra- 
fia e de jazigos minerais de elevado valor, nomeada- 
mente a última que se pode considerar excepcional. 

O Laboratório de Mineralurgia e Planeamento 
Mineiro (LMPMIST) está apetrechado para a reali- 
zação de ensaios laboratoriais descontinuos de mine- 
ralurgia cobrindo os processos fundamentais da Pre- 
“paração de Minérios. Dispõe ainda de certo equipa- 
mento semi-industrial que permite o funcionamento 
de diagramas parcelares em escala piloto, Quanto a 
cálculo automático, dispõe o LMPMIST de facilida- 
des avançadas na área da Geoestatística e Planea- 
mento Mineiro. 


HI 


O Curso de Engenharia de Minas suporta-se em 
uma base alargada à matemática, física, química e 
petrologia-geologia. É um curso cujas raizes têm de 
mergulhar, por um lado, nas fontes fisico-matemáticas 
comuns a todos os cursos de engenharia, e por outro 
têm de ser alimentadas por um profundo conheci- 
mento das ciências da natureza de base químico-mine- 
ralógica e geológica. 

É este um dos principais aliciantes do Curso, mas 
também uma das suas principais características por 
vezes de difícil apreensão pelos alunos recém-chegados 
ou para ele “empurrados” por um “numerus clausus” 
cego como o que vigora entre nós. 

O Curso de Minas do IST termina oferecendo duas 
opções de ramos: o de Geologia Aplicada e o do Pla- 
neamento Mineiro. O primeiro põe o acento tónico 
no campo da prospecção mineira e da geologia apli- 
cada. O segundo privilegia o planeamento mineiro e a 
mineralurgia. 


Quanto a cursos de pós-graduação está em funcio- 
namento o curso de Mestrado em Mineralurgia e 
Planeamento Mineiro. 

Impõe-se o lançamento de outros cursos de Mes- 
trado nomeadamente no campo da Prospecção mi- 
neira e da Geomecânica. Esperamos que em breve 
sejam realidade. 
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IV 


Já se referiu que junto ao LAMPIST funciona um 
centro de investigação suportado financeiramente pe- 
lo INIC — Centro de Petrologia e Geoquímica da 
UTL (CEPGUTL). 


Como nunca os laboratórios, museus e demais es- 
tabelecimentos anexos ao IST previstos em diploma 
de há 61 anos foram providos por verbas especificas 
do Orçamento Geral do Estado têm subsistido pelo 
dinamismo do corpo docente e de investigação neles 
trabalhando e por via de algumas receitas dos traba- 
lhos executados para fora do Instituto, como prevê o 
artigo 59 do Decreto n. 7727. 

Todavia tal situação não se compadece com a in- 
vestigação científica e tecnológica dos nossos dias, 
pelo que se propôs, oportunamente, ao INIC, a cria- 
ção do CEPGUTL. 

Fundamentalmente no Centro (e consequentemen- 
te no LAMPIST) hã duas grandes linhas de investi- 
gação. Uma ocupa-se da geoquimica da vulcanologia 
e da investigação geotérmica; outra da alteração e 
alterabilidade de rochas e das suas incidências geotéc- 
nicas. 

Na primeira linha de investigação tem-se procura- 
do fazer o estudo petrológico e geoquimico das ro- 
chas vulcânicas dos arquipélagos da Madeira, Selva- 
gens e das ilhas de S. Tomé e do Príncipe. Do mesmo 
modo se tem realizado investigação geotérmica, quer 
no domínio das altas entalpias (nos Açores), quer no 
das baixas entalpias (no continente). 

Os estudos realizados sobre as lavas madeirenses e 
das ilhas do Golfo da Guiné têm originado várias 
notas cientificas apresentadas nas diversas Semanas 
de Geoquímica que se têm realizado entre nós. Sobre 
o arquipélago das Selvagens, um docente do LAMPIST 
ultima o seu trabalho de tese de doutoramento. A 
colaboração com o Instituto de Geociências dos Aço- 
res, primeiro, e depois com o Laboratório de Geo- 
ciências e Tecnologia dos Açores, deu origem a vários 
estudos petrográficos, mineralógicos e geoquímicos 
relacionados com o Projecto Geotérmico dos Açores. 
Ainda no domínio da geotermia têm-se publicado no- 
tas sobre as potencialidades geotérmicas no domínio 
das baixas entalpias no continente. No recente Semi- 
nário sobre Tecnologias das Novas Energias (1982) 
realizado no âmbito do Estudo Tecnológico da In- 
dústria Portuguesa foi apresentado e discutido o rela- 
tório “Energia Geotérmica” da autoria de um investi- 
gador do CEPGUTL. Do mesmo modo foi o CEPGUTL, 
pelos seus investigadores, chamado a colaborar no 
Plano Energético Nacional (PEN), e no Grupo de 
Trabalho sobre | D&D sobre Energia. Ai foi dada 
particular atenção à participação da Energia Geo- 
térmica. 


Ainda no campo da colaboração com entidades 
exteriores ao LAMPIST/ CEPGUTL deve mencionar- 
-se a colaboração com os Serviços Geológicos de Por- 
tugal (ilha da Madeira, petrologia da folha de Lou- 
res), com a JICU/IICT (geoquímica das ilhas de S. 
Tomé e Príncipe), Instituto de Geociências da Uni- 
versidade Federal da Bahia (Salvador) em cursos de 
pós-graduação sobre geoquímica de vulcanitos, com 
a Ecole des Mines de Paris (petrologia da ilha da 
Madeira), com a Université de Paris-Sud (Dep. de 
Petrografia), nos estudos de rochas de S. Miguel 
(Açores), com a B.R.G.M. em especial com o Prof. 
H. de la Roche sobre cálculo geoquímico. No domi- 
nio das novas energias e da conservação da energia, 
investigadores do LAMPIST/CEPGUTL não só têm 
colaborado estreitamente em acções com o LGT dos 
Açores, como da Ordem dos Engenheiros, como ain- 
da em projectos realizados pela JNICT e muito espe- 
cialmente em recente “Workshop” com apoio da U.S. 
National Academy of Sciendes levado a cabo pór 
aquela Junta. Há ainda ligação institucional com a 
International Association of Volcanology and Che- 
mistry of the Earth's Interior. 

No que concerne à segunda linha de investigação 
— alteração e alterabilidade de rochas, suas incidên- 
cias na geotecnia — há bastante trabalho realizado. 
Este divide-se em trabalho experimental e em estudos 
de campo com apoio laboratorial. Do primeiro tipo 
de investigação resultou um lote de publicações em boa 
parte feitas em revistas internacionais da especialida- 
de (Engineering Geology e Bulletim da IAEG) ou 
apresentados nos Congressos Internacionais da So- 
ciedade Internacional de Mecânica das Rochas e da 
IAEG. Fundamentalmente tem-se procurado estudar 
o mecanismo da alteração e da alterabilidade de ro- 
chas (endógenas e exógenas), definindo-se índices de 
alterabilidade. Acresce que se tem procurado a apli- 
cação da investigação experimental laboratorial quer 
a problemas geotécnicos especificos, quer ao compor- 
tamento de rochas ornamentais, quer ainda no acompa- 
nhamento da análise dos fenómenos de degradação 
de rochas empregues em monumentos (Sé Velha de 
Coimbra e Mosteiro da Batalha). 

Há um docente trabalhando sobre a geoquimica da 
alteração dos perfis geológicos dos principais tipos 
litológicos da Serra de Sintra e outro dedicando-se ao 
estudo da alterabilidade dos calcários cristalinos da 
região de Borba-Estremoz-Vila Viçosa. Espera-se que 
ambos os trabalhos conduzam a teses de doutoramento. 

A colaboração com entidades exteriores ao IST no 
domínio desta linha de investigação tem sido estreita 
com a Direcção-Geral dos Edifícios e Monumentos 
Nacionais e com várias empresas exploradoras de 
mármores. Inclusivamente tem-se colaborado com a 
associação de industriais do ramo da exploração de 
pedras ornamentais (mármores e granitos) — 


ASSIMAGRA. No âmbito do Grupo de “Trabalho 
Português de Conservação de Pedras em Edifícios a 
colaboração tem sido intensa. A nível internacional 
tem-se colaborado com o Comité Petrographie do 
ICOMOS (UNESCO). 

Para além dos trabalhos inerentes às linhas de in- 
vestigação referidas colabora o LAMPIST/CEPGUTL 
em acções de serviços à comunidade (Petrografia, 
análises mineralógicas pelos raios-X, ATD, ATP, 
IV, preparação de amostras e minerografia de 
minerais opacos). Tem facultado estágios a investiga- 
dores de outros departamentos (Laboratório de Geo- 
ciências dos Açores, LNETI, Serviço de Estudos do 
Ambiente, IPT de S. Paulo, Brasil, Instituto de Geo- 
ciências da U.F. Bahia, Salvador, Instituto de Geo- 
ciências da U.F. do Rio de Janeiro). 

Dispõe o LAMPIST/CEPGUTL de uma razoável 
biblioteca provida com os principais livros e revistas 
do domínio da mineralogia, petrologia, e geoquímica. 
Esta biblioteca é complementada pela do Laborató- 
rio de Geologia nos campos da geologia, jazigos mi- 
nerais e geologia aplicada. 

Deve referir-se ainda a existência do Museu Ben- 
saúde. Trata-se de um museu de mineralogia e petro- 
logia especialmente portuguesas. A quantidade e a 
qualidade das espécies minerais e das amostras de 
rochas existentes é importante e tem possibilitado, 
desde a sua criação, estudos publicados relevantes. 

O pessoal científico trabalhando no CEPGUTL/ 
[| LAMPIST é de 10 investigadores, 6 técnicos de in- 
vestigação, preparadores e auxiliares de laboratório. 

O Laboratório de Geologia tem estado essencial- 
mente vocacionado para o ensino. No entanto, o seu 
corpo docente tem desenvolvido sempre actividades 
investigadoras embora não ligado a qualquer centro 
de investigação. 


Só para citar estudos recentes, assinala-se: 


— Investigação paleontológica (faunas paleozóicas). 

— Investigação geofísica (métodos geoeléctricos). Já deu 
lugar a uma dissertação para doutoramento. 

— Investigação geohidrológica (sistemas fracturados). 
Está em preparação uma dissertação para douto- 
ramento. 

— Investigação metalogenética (jazigos minerais me- 
tálicos). 


Toda esta investigação tem dado lugar a publicações. 


As nove disciplinas a cargo do Grupo de Explora- 
ção, cobrindo assuntos de Geomecânica, de Explora- 
ção de Minas e de Máquinas (especialmente minei- 
ras), não obstante o seu carácter essencialmente 
formativo, estão dirigidas para as aplicações em tra- 
balhos de campo e na indústria, a uma escala muito 
diferente da correspondente a uma investigação de 
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indole tecnológica que possa ser desenvolvida facil- 
mente na Escola. As características próprias daque- 
las matérias têm dificultado até a aquisição do equi- 
pamento laboratorial desejável, cujo elevadíssimo 
custo não se ajusta à satisfação de necessidades es- 
sencialmente pedagógicas; por esta razão, as aulas 
práticas laboratoriais de Geomecânica têm sido da- 
das no Laboratório Nacional de Engenharia Civil, 
beneficiando de Convénio com aquela instituição. 

O exposto não significa que, sempre que possível, O 
pessoal docente ligado ao Grupo de Exploração não 
exerça actividade de investigação, quer em colabora- 
ção com o Laboratório de Mineralogia e Petrologia 
do IST quer com outros organismos oficiais, tais co- 
mo o Laboratório Nacional de Engenharia Civile o 
Serviço de Fomento Mineiro. Dessa actividade são 
testemunhos numerosos relatórios referentes a traba- 
lhos específicos e um número também já elevado de 
publicações em Portugal e no estrangeiro, com ênfase 
para as comunicações apresentadas em reuniões cien- 
tíficas internacionais. 


No Laboratório de Mineralurgia e Planeamento 
Mineiro cursam-se disciplinas ao nível do 4º e 5.º 
anos dos cursos de Engenharia de Minas e Metalur- 
gia nas áreas da Mineralurgia e Planeamento Mineiro. 

As disciplinas de Mineralurgia têm forte compo- 
nente laboratorial, As de Planeamento Mineiro re- 
correm intensamente a cálculo computorizado. 


No ano escolar de 1982/83 iniciou-se no LMPMIST 
o Primeiro Curso de Mestrado em Mineralurgia e 
Planeamento Mineiro, no qual se contou e conta com 
colaboração de docentes da Universidade de Leeds, 
do Laboratório de Mineralogia e Petrologia do IST, 
da Faculdade de Ciências e Tecnologia da Universi- 
dade de Coimbra € de especialistas diversos. 


Espera-se que os alunos que frequentam este curso 
de mestrado venham a apresentar as dissertações no 
ano escolar de 1984/85. 


No ano escolar de 1981/82 o Laboratório de Mine- 
ralurgia e Planeamento Mineiro ministrou curso de 
técnicas laboratoriais a funcionários da Direcção Na- 
cional de Geologia da R.P. de Moçambique. 

Quer na área da Mineralurgia, quer na do Planea- 
mento Mineiro, o Laboratório tem dado colaboração 
a projectos mineiros nacionais. A contribuição pres- 
tada cobre principalmente domínios de conhecimento 
específicos não disponíveis noutras instituições esta- 
tais ou privadas do Pais. 

Junto do Laboratório de Mineralurgia e Planea- 
mento Mineiro funciona o Centro de Valorização de 
Recursos Minerais (CVRM) do INIC., 

Duas linhas de investigação do CVRM desenvolvem- 
-se no Laboratório de Mineralurgia e Planeamento 
Mineiro: 
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Linha | — Mineralurgia 
Linha 2 — Planeamento Mineiro 


Na linha | tem sido dado realce à modelagem de 
processos mineralúrgicos, tema central de dissertação 
para doutoramento e, bem assim, à sua automação. 
Este último tema tem sido trabalhado por bolseiro do 
CVRM no estrangeiro. 

Na Linha 2 desenvolveu-se ultimamente o Projecto 
SAVARN subsidiado parcialmente pela JNICT e em 
conexão com a Universidade de Leeds. 

Consequentemente criou-se um sistema integrado de 
programas que permite o tratamento informático dos 
dados resultantes de campanhas de prospecção de re- 
cursos minerais (ou mesmo de outros recursos natu- 
rais) com vista à sua avaliação e valorização. 

Recentemente doutoraram-se dois docentes do 
LMPMIST na área do Planeamento Mineiro em co- 
nexão com a actividade desenvolvida pelo citado 
projecto. 

Foi também possivel instalar meios informáticos 
no LMPMIST que permitem o adequado tratamento 
de problemas do Planeamento Mineiro e está pro- 
posto novo projecto onde se pretende desenvolver a 
informática da simulação de processos mineralúrgi- 
cos e sua automação. 

Outro docente prepara o seu doutoramento apoia- 
do na investigação desenvolvida na Linha 2, no âmbi- 
to da Análise de Dados. 

O CVRM tem facultado numerosos contactos com 
centros de investigação e universidades estrangeiras. 
A contribuição dos docentes do LMPMIST em con- 
gressos e reuniões similares e em revistas de especiali- 
dade, quer nacionais, quer estrangeiras, tem sido 
muito satisfatório nos últimos dez anos. 


V 


Os docentes do Curso de Minas, de todos os esca- 
lões, têm colaborado intensamente com a Indústria, 
com Departamentos de Investigação Aplicada da 
Administração pública, com a Ordem dos Engenhei- 
ros, etc.. Acresce que participam nas mais variadas 
reuniões técnico-científicas nacionais ou internacio- 
nais como seminários, congressos, cursos avançados, 
cic.. 

Merece menção especial a colaboração estreita 
com as principais companhias mineiras, inclusivamen- 
te realizando os mais variados tipos de trabalhos, 
desde o estudo mineropetrográfico de sondagens até 
à participação no planeamento de métodos de lavra e 
ensaios semi-industriais em lavaria-piloto. 

A colaboração com laboratórios nacionais como o 
LNEC (no campo da geomecânica), o LNETI (no 
domínio da mineralurgia e da energia) e com o Insti- 


tuto de Investigação Científica Tropical (IICT / JICU) 
tem sido relevante. 

Em relação à Ordem dos Engenheiros, a Especiali- 
dade de Engenharia de Minas é fortemente suportada 
pela actividade dos docentes do IST. Inclusivamente 
a organização anual dos Encontros Nacionais de En- 
genharia de Minas (de que se avizinha a realização do 
4.º, conta sempre com o impulso decisivo dos docen- 
tes e antigos alunos do Curso de Minas do IST que 
confraternizam com os colegas provenientes de outras 
escolas superiores de engenharia. Em outras manifes- 
tações orientadas pela Ordem dos Engenheiros se tem 
notado a participação efectiva da nossa Escola em 
geral e dos docentes de Minas em especial (Jornadas 
de reflexão sobre o Plano Energético Nacional, Jor- 
nadas de engenharia dos países lusófonos, em organi- 
zação, Congressos da Ordem, etc.). 

Acresce que os docentes de minas têm sido chama- 
dos a dar colaboração especializada, técnico-cientifica, 
a vários Conselhos Superiores e Consultivos. Assim 
têm assento no Conselho Superior de Engenharia e 
Tecnologia Industrial, no Conselho Superior de Mi- 
nas, no Conselho Superior de Obras Públicas, no 
Conselho Científico do Instituto Nacional de Investi- 
gação Cientifica, etc. 

Na nossa Escola e nos vários Laboratórios do Cur- 
so de Minas se têm realizado algumas manifestações 
científicas de que se referem apenas a 1.º Semana de 
Geoquímica promovida pelo LAMPIST em 1977 e 
que representou o lançamento destas reuniões no 
país. O MPMIST promoveu o 1.º Seminário de 
Geoestatiística em 1978 consistindo também um mar- 
co indelével na investigação geomatemática no pais. 

Nas suas relações com laboratórios congêneres es- 
trangeiros tem havido amplo intercâmbio de investi- 
gadores. São de realçar as relações com as Universi- 
dades brasileiras (UF Bahia, Rio de Janeiro e U. de S. 
Paulo e ainda com o IPT de S. Paulo), com as Uni- 
versidades de Luanda e Eduardo Mondlane (Mapu- 
to), e com universidades e demais departamentos 
europeus e norte-americanos (École des Mines de Pa- 


ris, Bureau de Recherche Géologique et Miniêres, 
Universidade de Minesota, Cambridge, Leeds, Li- 
bre de Bruxelles, de Pisa. 


VI 


Procurando dar uma visão prospectiva de domi- 
nios dos recursos naturais minerais e energéticos no 
pais que parecem promissores para actividades de 
ID&D, apresenta-se a seguinte lista genérica: 


[D) Cartografia geo-mineira e prospecção geoqui- 
mica e geofísica. 
ll) Cartografia geotécnica e análise geomecânica 
das fundações de grandes obras de arte e de zonas 
a urbanizar, de parques industriais, etc.. 

WI) Cartografia geohidrológica de áreas sensíveis do 
país. Avaliação da capacidade dos recursos geo- 
hidrológicos e sua exploração racional para fins 
múltiplos. 

IV) Metalogenia do cobre, volfrâmio e estanho. 

V) Geomecânica aplicada à lavra mineira e tecno- 
logia de sondagens profundas. 

VI) Exploração racional de rochas, quer como ma- 
teriais de construção, quer como matérias-primas 
(areias, argilas, diatomitos) quer como rochas or- 
namentais. 

VII) Emprego racional e conservação de pedras em 
infra-estruturas (portos, vias de comunicação, 
etc.) e em edifícios e monumentos. Alteração e 
alterabilidade de rochas. 
VIII) Planeamento mineiro (v.g. grandes projectos co- 
mo o do ferro, do carvão, das pirites, etc.). 
IX) Tratamento mineralúrgico. 
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MAPAS DE MECANISMOS COMPETITIVOS 
DE RESPOSTA DUM SISTEMA 


M. AMARAL FORTES Departamento de Metalurgia, Instituto Superior Técnico. / 
Centro de Mecânica e Materiais da U.T.L. (CEMUL) 


RESUMO 


Descrevem-se sumariamente os mapas de meca- 
nismos já utilizados em Ciência dos Materiais e 
indicam-se as suas vantagens. Estes mapas são de 
aplicabilidade mais geral; discutem-se as situações em 
que podem ser construídos e faz-se a sua classificação 
consoante a natureza das variáveis que intervêm. Fi- 
nalmente, dão-se exemplos de mapas de mecanismos 
que poderiam ser elaborados para descrever fenóme- 
nos e processos fora da Ciência dos Materiais. 


ABSTRACT 


Mechanism maps in use in Materials Science are 
described and their advantages are outlined. Maps of 
this type are of general applicability; we discuss the 
situations in which they can be constructed and the 
various categories in which they may fall. Various 
examples of maps outside Material Science are 
presented. 


1. MAPAS DE MECANISMOS DE DEFORMAÇÃO 


Em 1972 foi publicado na revista Acta Metallurgica 
um artigo de M.F. Ashby (Engineering Department, 
Universidade de Cambridge) com o título “A first 
report on deformation mechanism maps” (1). Nele se 
apresentava um método extremamente simples e ino- 
valor de descrição dum fenómeno muito complexo 
— a deformação plástica de sólidos. A complexidade 
da deformação dos metais, dos cerâmicos e vidros e 
dos polimeros, resulta da grande variedade de meca- 
nismos que contribuem para a deformação e da difi- 
culdade em definir a relevância relativa desses meca- 
nismos. Por exemplo, no caso dos materiais metálicos 
e cerâmicos, policristalinos (englobando os aços e as 
ligas de alumínio, o gelo e os óxidos metálicos refrac- 
tários), é possível identificar cerca de dez mecanismos 
de deformação, os quais têm a ver com o tipo de 
defeitos móveis responsáveis pela deformação (deslo- 
cações, defeitos pontuais) e com a forma como eles se 
movem sob a acção de tensões aplicadas (por onde se 
movem e quais os obstáculos com que deparam). 
Conforme o valor da temperatura, T, e da tensão 
aplicada, o (tracção duma barra cilíndrica, por 
exemplo) assim prepondera um ou outro mecanismo, 

Os mapas de mecanismos de deformação são gráfi- 
cos, com eixos (T,o ), onde estão assinaladas as re- 
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giões em que cada mecanismo é dominante. Os ma- 
pas contêm também informação sobre a velocidade 
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FIGURA 1 — Mapa de deformação com eixos tensão, o, € 
temperatura, T, para o níquel de tamanho 
de grão d =0,1] mm. As linhas grossas mar- 
cam as fronteiras das regiões em que pre- 
domina cada um dos mecanismos indica- 
dos. As linhas finas são linhas de velocidade 
de extensão (É) constante. G é o módulo 
de rigidez e T; a temperatura absoluta de 
fusão. 
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FIGURA 2 — Mapa de deformação o — T para o gelo de 
tamanho de grão d= 10 mm. Ver legenda 
da Fig. 1. 
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família de curvas é = constante traçadas no mapa. 
Nas Figs. | e 2 dão-se exemplos de mapas de defor- 
mação, respectivamente para o níquel e para o gelo (2). 
Apenas um exemplo: a 400º C e sob uma tensão de 
102 MN / m'o níquel deforma-se com uma velocidade 

= 10 *s”! e o mecanismo preponderante de de- 
formação é a fluência de baixa temperatura, envol- 
vendo o deslizamento e trepa de deslocações. 

Quaisquer duas das três variáveis independentes 
o, T, é, podem ser usadas para eixos do mapa. O 
mapa da Fig. 3 contém a mesma informação que o da 
Fig. 1, mas as coordenadas do mapa são agora €, T. 
Na maior parte dos mapas de deformação, usa-se o, | 
como coordenadas. 

São várias as vantagens destes mapas, autênticos 
catálogos do comportamento mecânico dos mate- 
riais. Para além da informação que contêm, uma das 
principais vantagens é, sem dúvida, o seu valor didác- 
tico. Com eles torna-se simples compreender a com- 
petição entre os diversos mecanismos de resposta 
dum material a esforços aplicados. 

A construção de um mapa de deformação para um 
dado material não é, no entanto, fácil. São necessá- 
rios os seguintes passos; 


| — Identificar os diversos mecanismos de defor- 
mação. 
2 — Obter, quer por via teórica quer empiricamente, 


as equações que descrevem o contributo de ca- 
da mecanismo para a deformação. Estas equa- 
ções relacionam a velocidade de extensão € 

com a tensão aplicada o ea temperatura I,e 
fazem intervir variáveis definidoras do estado 
“mecânico” do material (p. ex. a densidade dos 
defeitos responsáveis pela deformação e a den- 
sidade de obstáculos ao seu movimento) e ainda 
diversas propriedades do material, como módu- 
los de elasticidade, coeficientes de difusão, etc. 
Este é o passo dificil. 


3 — Uma vez de posse das equações descritivas de 
cada mecanismo, “basta” calcular a velocidade 
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FIGURA 3 — Mapa de deformação com eixos é — T para 
o níquel de tamanho de grão d =0,] mm, 
As linhas finas são linhas de tensão norma- 
lizada (o/G) constante. Comparar com o 
mapa da Fig. 1. 


de extensão é para cada um deles e ver, para 
cada o, T, qualé o mecanismo mais rápido. É 
esse o mecanismo dominante e, em primeira 
aproximação, o valor observado para é é ova- 
lor devido a esse mecanismo, 


É preciso ainda fazer algumas hipóteses simplifica- 
tivas para poder chegar ao fim. A principal é a que 
consiste em admitir que o estado mecânico do mate- 
rial não se altera por deformação ou, o que é menos 
restrictivo, que em cada regime se atinge um estado 
estacionário. A inclusão de variações de estado mecã- 
nico e de transientes pré-estacionários complica a 
elaboração dos mapas e obriga a introduzir outras 
variáveis (2); de resto faltam, em regra, as equações 
descritivas dessas situações, 

O êxito dos mapas de deformação tem sido muito 
grande. Até agora foram publicados mapas para cer- 
ca de 50 materiais metálicos'e cerâmicos. O mais re- 
cente é o mapa de deformação do cobalto (3). Uma 
excelente compilação destes mapas apareceu em 1983 
(ref. 2). 


2. OUTROS MAPAS EM CIÊNCIA 
DOS MATERIAIS 


Como é natural, tendo em vista o êxito dos mapas 
de deformação, os mapas de mecanismos alastraram 
a outros sectores da Ciência dos Materiais com o 
intuito de descrever outros fenómenos complexos em 
que intervêm diversos mecanismos. Foram propostos 
e elaborados mapas de sinterização (4,5) e mapas de 
fractura (6). Em ambos os casos os mapas mostram 
os dominios das variáveis em que prevalece este ou 
aquele mecanismo e dão informação sobre a cinética 
global desses fenómenos. 

Nos mapas de sinterização (Fig. 4) as variáveis são 
a densidade e a temperatura. A cinética é indicada 
por linhas At =constante, sendo Ato tempo reque- 
rido para se atingir uma determinada densidade a 
uma determinada temperatura, a partir dum estado 
em que não ocorreu ainda sinterização do pó. Repare-se 
que uma das variáveis (a densidade) é uma proprie- 
dade do sistema e não uma variável que se possa 
“controlar exteriormente”. 

Os mapas de fractura deparam com diversas difi- 
culdades, que vão desde a dificuldade de escolha das 
variáveis e com a dificuldade de individualização dos 
mecanismos e passam pela falta de equações adequa- 
das a cada mecanismo. Estas dificuldades estão ob- 
viamente relacionadas com a nossa deficiente com- 
preensão do fenómeno da fractura. A elaboração de 
mapas para um dado fenómeno põe a nu as limita- 
ções do conhecimento sobre esse fenómeno — outra 
vantagem dos mapas. 
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Nas secções seguintes procuraremos definir as re- 
gras e os requisitos que estão por detrás da elabora- 
ção de mapas do tipo dos que foram descritos ante- 
riormente. É evidente que a sua aplicabilidade não se 
restringe à Ciência dos Materiais e procuraremos dar 
exemplos em outros domínios, de situações e fenó- 
menos susceptíveis de representação em mapas. Gos- 
tariamos muito de ter conhecimento de outros exem- 
plos, provavelmente melhores, que porventura ocor- 
ram aos leitores deste artigo. De resto, o autor não 
põe de lado a hipótese de mapas deste tipo serem já 
usados em outros domínios para além da Ciência dos 
Materiais. 


3. REGRAS E REQUISITOS PARA 
ELABORAÇÃO DE MAPAS 


Em termos gerais, a construção de um mapa de 
mecanismos é, em princípio, possível (e será portanto 
desejável) em todas as situações em que haja vários 
mecanismos alternativos ou paralelos conducentes ao 
mesmo resultado. 

Tentemos agora precisar os termos deste enuncia- 
do, inspirando-nos nos mapas descritos nas secções 
anteriores. Em todos os casos, há um sistema físico, 
macroscópio, que evolui sob a acção de um ou mais 
estímulos. Um estimulo pode traduzir um modo de 
interacção do sistema com o exterior (p. ex. campos 
aplicados) ou um mero desvio duma variável de esta- 
do do sistema (p. ex. a temperatura) em relação ao 
valor inicial de equilíbrio. Os estímulos podem ter 
diversos níveis de intensidade. Em consequência dos 
estímulos verifica-se um “acontecimento” no sistema, 
a que chamaremos resposta do sistema. À resposta 
pode ter vários níveis de intensidade ou vários niveis 
de rapidez (medida por ex. pelo tempo ao fim do qual 
ocorre um acontecimento singular). Finalmente ad- 
mitiremos que a evolução do sistema, e portanto a 
resposta, pode ocorrer por diversos processos ou me- 
canismos (envolvendo os átomos, moléculas, 10es, 
electrões, etc. do sistema) que têm por palco o siste- 
ma. Estes mecanismos podem ser de duas categorias: 
a) alternativos — só um ocorre em cada regime; b) 
paralelos — podem ocorrer simultaneamente mas a 
cinética de cada um deles é independente da existên- 
cia dos outros. Uma mistura de dois mecanismos À e 


FIGURA 4 — Mapa de sinterização de partículas esféricas 
de prata de 107” m de diâmetro. As coor- 
denadas são a densidade relativa do pó, 
ô/ôp, € a temperatura. As linhas mais 
grossas delimitam as regiões em que cada 
mecanismo de difusão (pelo interior das 
partículas, etc.) é preponderante na den- 
sificação. As linhas tracejadas marcam di- 
versas fases durante a sinterização. Por ex- 
emplo, na fase |, as partículas são ainda 
identificáveis e na fase 3 os poros estão Iso- 
lados. Estão traçadas também linhas À t= 
constante (ver texto). 


B (por exemplo, A e B em série) deve portanto ser 
considerada como um mecanismo adicional, €C. A 
temperatura tem o carácter de um estímulo, alteran- 
do a cinética dos diversos mecanismos, mas convém 
considerá-la explicitamente. 

Há vantagem em distinguir dois casos no que diz 
respeito à evolução do sistema sob a acção dos esti- 
mulos e à forma como se quantifica a resposta: 


| — Osistema evoluie a resposta é caracterizada pelo 
tempo para atingir um valor especificado duma 
propriedade; ou ainda pelo tempo que decorre 
até que se verifique um acontecimento singular 
ou que se atinja um estado de equilíbrio, 


2 — Como caso especial do anterior, o sistema poderá 
atingir sob a acção dos estimulos um valor esta- 
cionário duma propriedade, valor esse que po- 
de ser utilizado para medir a resposta. 


Poder-se-ia ainda admitir que a resposta assume 
aspectos qualitativamente diferentes, conforme o me- 
canismo preponderante (ver exemplo na secção 6.2). 

Finalmente, note-se que o termo “propriedade” de- 
ve ser tomado num sentido lato, de modo a incluir a 
taxa de variação duma propriedade, ou um fluxo, ou 
qualquer grandeza característica dum estado estacio- 
nário ou não. 

Como se pretende que os mapas sejam bidimensio- 
nais, há que restringir o número de variáveis. Se se 
pretender representar o efeito de variáveis internas do 
sistema (p. ex. o tamanho de grão em sólidos policris- 
talinos, no caso da deformação ou fractura) poderá 
não haver outro remédio senão elaborar mapas para 
diversos valores dessas variáveis (mapas tridimensio- 
nais são de dificil construção e leitura). 

Sintetizando: se for E a variável que mede a inten- 
sidade do estimulo (consideraremos a existência de 
apenas um estimulo), T a temperatura, V, as variá- 
veis de estado com relevância nos diversos mecanis- 
mos e Pg os parâmetros característicos do sistema, 
independentemente do estado em que se encontra (p. ex, 
a estrutura cristalina dum sólido, energias de activa- 
ção, etc.), a resposta R! para cada mecanismo isola- 
do é função de todas estas variáveis: 


Ri=f; (E, T; Va: Pp) [1] 


Esta função tem de ser obtida para cada mecanismo 
quer por via teórica quer empiricamente. Neste últi- 
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mo caso é preciso averiguar se a equação se pode 
extrapolar para a gama de valores das variáveis em 
que outros mecanismos predominam. A resposta Ri 
pode ser um intervalo de tempo (At),, ou uma 
grandeza com o carácter dum fluxo ou de taxa de 


. Mr E. E = r F 
variação duma caracteristica do sistema ts 
t 


Se se conhecerem as equações do tipo [1] para 
cada mecanismo é possivel calcular, para cada con- 
junto de valores das variáveis, qual é o mecanismo 
dominante. Como se supõe que os mecanismos são 
paralelos, o mecanismo dominante é sempre o mais 


rápido (menor (At), ou maior valor absoluto des ). 


Para construir mapas bidimensionais é preciso se- 
leccionar duas variáveis apenas (além da resposta) e 
considerar que as restantes se mantêm constantes. Só 
pode haver portanto um estímulo variável além da 
temperatura; a terceira variável caracteriza a respos- 
ta. Se o único estimulo for a temperatura, a segunda 
variável será uma das variáveis internas Vá. por 
exemplo V,. A consideração dos parâmetros Pg só 
interessa quando se pretendem comparar classes de 
sistemas (p.ex., o comportamento de diversos metais). 

As situações possiveis caem portanto num dos se- 
guintes tipos. 


Tipo 1) RE=HÃ(E, T) Va = constante 

em que os eixos no mapa são E e T. É útil traçar no 
mapa linhas R =const. em que R é a resposta global, 
observável. Claro que qualquer outro par de variáveis 
(E, R)ou(R, T) podem ser usadas como coordenadas 
do mapa. 


Tipo 1) Ri=fi(T.V)) Va= cont a > 1 

em que os eixos do mapa são, por exemplo, Vje T: 
é útil traçar no mapa linhas R = const. À varia- 
ção de V, pode ocorrer durante a evolução do sis- 
tema, ou pelo simples facto de o seu valor inicial 
poder ser variado. Se todas as resposta R! forem 
proporcionais ao estímulo E, podem fazer-se mapas 
do tipo II considerando sempre um estímulo unitário. 

No caso da resposta poder ser qualitativamente di- 
ferente conforme o valor das variáveis, torna-se ne- 
cessário arranjar um critério que permita definir o 
nível de resposta, independentemente do tipo de 
resposta, 

Resta a seguinte questão: como se obtem o valor R 
da resposta global a partir dos valores Ri das respos- 
tas correspondentes a cada mecanismo, supondo-o 
isolado? Se não for conhecida a lei de composição 
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R=f(Ri) 
pode admitir-se a aditividade das respostas: 


R=> Ri 


ou, no extremo da simplificação, admitir-se que a 
resposta global tem o valor da resposta correspon- 
dente ao mecanismo dominante (R'q 


R=(Ri')g 


4. RESPOSTA CONTROLADA POR 
OBSTÁCULOS 


Um caso particular do tratamento anterior que 
convém ser explicitado, pela frequência com que se 
depara, é aquele em que a resposta é dificultada por 
vários obstáculos de natureza diferente, em série. O 
poder impeditivo dos obstáculos varia com o estimu- 
lo e com a temperatura. É o que acontece com os 
obstáculos ao movimento de transportadores de car- 
ga (dum só tipo): a barreira de energia tem uma altu- 
ra efectiva que depende do estimulo (um campo eléc- 
trico aplicado) e a passagem da barreira é feita por 
activação térmica. Para densidades fixas de obstácu- 
los (os Va da eg. 1), pode traçar-se um mapa em que 
as coordenadas são E e T e em que se assinalam as 
regiões em que cada obstáculo é mais “poderoso”. 
Podem traçar-se linhas de velocidade média constan- 
te dos transportadores e admitir que há aditividade nos 
tempos de passagem dos vários tipos de obstáculos. 


5. IDENTIFICAÇÃO DOS MAPAS USADOS 
EM CIÊNCIA DOS MATERIAIS 


Consideremos, de novo, os mapas referidos nas 
secções le 2 a fim de exemplificar o formalismo geral 
desenvolvido na secção 3. 


5.1 — Mapas de Mecanismos de Deformação 


O estímulo E é a tensão aplicada o. A resposta é 
medida pela velocidade de extensão, é. As coorde- 
nadas do mapa são T, o etraçam-se linhas é = const. 
Admite-se que o sistema atinge um estado estacioná- 
rio em cada regime, variável de regime para regime, 
em que se mantem constante a densidade de defeitos 
móveis. Esta hipótese é razoável às temperaturas 
mais altas. Para menores temperaturas não há esta- 
cionaridade e supõe-se que o valor de é é o corres- 
pondente à resposta inicial. Estes mapas são do tipo 1; 
o mecanismo dominante é o que origina maior e€. 

O tamanho de grão é um exemplo de variável que 


interfere na forma dos mapas e que se supõe ter um 
valor fixo em cada mapa. Se se quiser usar o tamanho 
de grão como variável do mapa, há que admitir que a 
tensão (p.ex.) tem um valor fixo (8). Neste caso o 
mapa será do tipo II, muito embora se admita que o 
tamanho do grão não varia com a deformação. 


5.2 — Mapas de Fractura 


O estimulo pode ser a tensão aplicada (6) ou outra 
grandeza mais adequada que tenha em conta o tama- 
nho inicial da fissura macroscópica responsável pela 
fractura, p.ex., o chamado factor intensidade de ten- 
são (7). As outras variáveis são a temperatura e o 


tempo para fractura ou a velocidade e de aumento 
de comprimento, a, da fissura. Traçam-se linhas de tem- 


E da 
po para fractura constante ou de dç constante. Estes 


mapas são do tipo I: o mecanismo controlador é o que 


origina menor (At); ou maior ca ; 


5.3 — Mapas de Sinterização 


Estes mapas exemplificam o tipo II. O estado ini- 
cial do pó não sinterizado é um estado de não equili- 
brio. A eliminação de porosidade ocorre por vários 
mecanismos que dependem de T e da densidade ô em 
cada instante. O mecanismo controlador é o que ori- 


gina o maior valor de Ea ou seja da taxa de densi- 
ficação. As variáveis do mapa são de T e traçam-se 


linhas ce = const. ou, se se preferir, linhas A t= const. 
sendo At o tempo necessário para atingir (à tempe- 
ratura T) a densidade ô a partir do estado inicial não- 
-sinterizado. 


6. OUTROS EXEMPLOS DE MAPAS 


Nesta secção apresentam-se exemplos de outros fe- 
nómenos e processos susceptíveis de representação 
sob a forma de mapas de mecanismos. Os mapas se- 
rão descritos duma forma esquemática sem pretensão 
de os construir para sistemas reais. 


6.1 — Mapas de Mecanismos de Transporte 


Suponhamos que há vários tipos de transportado- 
res de carga (ou de calor, ou de matéria). Seja p; a 
densidade, p; o valor da carga (ou da propriedade 
que for transportada) dos transportadores ie M; a 


respectiva mobilidade. O estímulo E poderá ser a in- 
tensidade dum campo eléctrico (ou um gradiente de 
temperatura ou de concentração). À corrente |; trans- 
portada pelos transportadores 1 é definida como a 
quantidade de propriedade que passa por unidade de 
tempo numa unidade de área, e é dada por 


=p; PivVi 


em que v; é a velocidade dos transportadores 1, rela- 
cionada com a mobilidade M; por 


vi= Mi; (E, 1). E 


Em geral M; é função de E e de T. E em geral p; é 
função de T podendo igualmente ser função de E. As 
equações descritivas são do tipo 


Il; = fi (E, T) [2] 


e a resposta global é 
| = > l; 
l 


Conhecida a forma das eqgs. [2] é possível assinalar 
nos mapas (coordenadas E, T) as regiões onde pre- 
domina cada um dos transportadores e traçar linhas 
|=constante. O mapa é do tipo 1. 


6.2 — Mapas de Transformação de Fase 


Suponhamos que uma fase, estável numa dada 
gama de temperatura acima de T. é arrefecida brus- 
camente até uma temperatura T < T.. Haverá uma 
transformação de fase, com uma cinética que depen- 
de de T e originando novas fases que dependem tam- 
bém de T. O estímulo para a transformação é a va- 
riação de energia livre associadaa AT=T.-T. Se 
todas as outras variáveis forem fixadas é possível tra- 
çar mapas com coordenadas T (ou AT) et, o tempo 
para que se transforme uma dada percentagem da 
fase inicial. As curvas respectivas designam-se por 
curvas [TT (IT) e são muito usadas para programar 
os tratamentos térmicos dos aços. Trata-se de um 
mapa de mecanismos em que a resposta é qualitati- 
vamente diferente conforme o valor do estímulo. O 
nível de resposta pode ser indicado pela percentagem 
da fase inicial transformada. Estes mapas não são, 
porém, do tipo dos considerados até agora. Falta 
uma variável. Uma variável adicional que poderá ser 
importante considerar é qualquer das seguintes: a) o 
tamanho de grão da fase inicial, d; b) o valor duma 
tensão mecânica, o, aplicada durante a transforma- 
ção à temperatura T; c) a velocidade de arrefecimento 


sl desde a temperatura inicial até T. Os mapas con- 
t 
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terão como variáveis, além de T e t (tempo para que 
se decomponha uma certa fracção da fase inicial) a 
terceira variável a), b) ou c). Em cada região indica-se 
o produto de decomposição da fase inicial. 


6.3 — Mapas de Mecanismos duma Reacção Química 


Consideremos uma reacção quimica que conduz à 
formação de um produto de reacção com um rendi- 
mento que é influenciado pela temperatura T e pelo 


pH do meio (p.ex.). Há vários mecanismos possíveis, 


envolvendo por exemplo estados intermédios diferen- 
tes; cada um deles tem uma cinética que é influencia- 
da por T e pelo pH. O mapa respectivo terá estas 
coordenadas e indicará os dominios em que ocorre 
um ou outro dos mecanismos. Podem traçar-se linhas 
de rendimento constante. 


6.4 — Mapas de Mecanismos de Morte 


Com este exemplo e com o seguinte pretende-se 
mostrar a aplicabilidade dos mapas a situações que 
não envolvem um sistema físico. Suponhamos que há 
uma população de bactérias num meio à temperatura 
Te com uma caracteristica quimica que pode ser va- 
riada (p.ex. a concentração, c, duma substância tóxi- 
ca). As bactérias podem morrer por uma de várias 
causas (mecanismos). A taxa de mortalidade para ca- 
da uma das causas depende de T, c, que serão as 
variáveis do mapa. Assinalam-se as regiões onde cada 
causa de morte é preponderante e traçam-se curvas de 
tempo médio de vida, por exemplo. É um mapa do 
tipo II. 


6.5 — Selecção de Processos 


Suponhamos que em vez de mecanismos diferentes 
há a considerar, mais propriamente, processos dife- 
rentes de obter o mesmo objectivo. Como exemplo, 
consideremos a extracção dum metal a partir de mi- 
nérios com caracteristicas diferentes, podendo recorrer- 
-se a vários processos de extracção e purificação. 

Pretende-se obter o metal com um grau de pureza 


Ta 


global determinado e com um teor duma determina- 
da impureza I não superior a um certo valor. Os mi- 
nérios de partida têm concentrações diferentes do me- 
tal, cm, e da impureza indesejável, cj. São estas as 
variáveis do mapa. O custo (por unidade de metal 
produzido) de cada processo de extracção é função de 
cm e cy. e haverá vantagem, por razões de preço, 
em usar um ou outro processo. Esta indicação é dada 
no mapa onde se podem traçar também linhas de 
preço unitário do metal extraido pelo processo mais 
económico. 


7. COMENTÁRIOS FINAIS 


A possibilidade de construção de “bons” mapas de 
mecanismos para um dado fenómeno, num dado sis- 
tema, implica um conhecimento teórico muito com- 
pleto do fenômeno e o conhecimento, por via expe- 
rimental ou teórica, de um certo número de caracteris- 
ticas do sistema que intervêm no fenómeno. À precisão 
ou qualidade dos mapas é portanto um bom aferidor 
do estado de conhecimento sobre o fenómeno em 
causa. Independentemente do seu poder de sistemati- 
zação e de aferição, os mapas, mesmo que pouco 
precisos por falta de dados suficientes, têm um valor 
didáctico notável. A sua aplicabilidade parece ser 
muito geral, mesmo fora das ciências físicas e da en- 
genharia, como se pretendeu mostrar neste artigo. 


REFERÊNCIAS 


(lj — M.F, Ashby, Acta Metall. 20, 887 (1972). 


(2) — HJ. Frostand M.F. Ashby,“Deformation Mechanism Maps”, 
Pergamon Press, 1983. 

(3) — P,M. Sargent, G. Malakondaiah and M.F. Ashby, Scripta Me- 
tall. 17, 625 (1983). 

(4) — M.F. Ashby, Acta Metall. 22, 275 (1974). 

(5) — F.B, Swinkels é M,F. Ashby, Acta Metall. 29, 259 (1981). 

(6) M.F. Ashby, C. Gandhi and D.M.R. Taplin, Acta Metall, 27, 
699 (1979). 

(7) — D,G. Brandon, Scripta Metall, 14, 555 (1980). 

(8) — T.G. Langdon and F.A, Mohamed, Materials Sci. Eng. 32, 
103 (1978). 


A EQUAÇÃO RELATIVISTA DO MOVIMENTO 
DE UMA BARRA NÃO RÍGIDA 


A. BROTAS E J.C. FERNANDES 
Instituto Superior Técnico 


RESUMO: 


As fórmulas da elasticidade relativista dos sólidos 
não rigidos apresentadas num anterior artigo (Técni- 
can: 461) são utilizadas no estabelecimento da equa- 
ção do movimento de uma barra. A passagem ao 
limite c— e conduz à equação clássica. 

É comentada a quase total ausência de equações do 
movimento dos sólidos deformáveis nos livros de 
Relatividade. 


SUMMARY: 


The formulae of relativistic elasticity of non-rigid 
presented in a previous paper (Técnica n. 461) are 
here used to derive the movement equation of a non- 
-rigid bar. When the limit c — ee is taken the classical 
equation in obtained. 

We comment on the almost omission in relativity 
books of equations dealing with the movement of 
deformable solids. 


Embora seja de registar o aparecimento logo em 
1911 de um importante artigo sobre a Elasticidade 
Relativista [ 1], esta matéria foi, sem dúvida, nos 70 
anos que se seguiram, o capitulo pobre dos estudos 
relativistas. 

A ausência de problemas científicos ou técnicos 
que se imediato exigissem a consideração de corpos 
sólidos relativistas (entenda-se modelos relativistas de 
corpos sólidos) e dificuldades conceptuais relaciona- 
das com estes corpos (que muito autores consideram 
inultrapassáveis no campo da Relatividade Restrita) 
fizeram com que a elasticidade dos sólidos fosse es- 
quecida (ostensivamente ignorada) na quase totali- 
dade dos livros e tratados de Relatividade. 

Surgiram, porém, nos últimos anos, indícios de 
uma mudança nesta atitude generalizada de voltar as 
costas aos corpos sólidos. 

Há, em particular, a assinalar a publicação, em 
1978, de dois livros relacionados com o assunto: “Re- 
lativistic Theories of Materials” de Aldo Bressan [2]e 

“Special Relativity” de W.G. Dixon [3). 

O livro de Bressan (que nos chegou às mãos há 
cerca de dois meses) é, segundo cremos, o primeiro e 
possivelmente até hoje único livro de Relatividade 
integralmente dedicado ao estudo dos meios conti- 
nuos “materiais” (1) sem exclusão dos sólidos (2). 


Original recebido para publicação em 13/1/81 


Universidade Técnica de Lisboa 


No livro de Dixon (em que os corpos sólidos ainda 
quase não aparecem) há a opção de tratar os pro- 
blemas dos meios contínuos só no âmbito da Relati- 
vidade Restrita, evitando a fuga para diante da ida 
para a Relatividade Generalizada, característica de 
tantos textos relativistas, que, a nosso ver, (neste ca- 
so) só complica os textos e nos afasta do problema de 
fundo. 


A procura de convenientes relações constitutivas 
que permitam determinar no referencial próprio as 
componentes do tensor impulsão-energia Te 8 dos 
meios contínuos considerados é, com efeito, a ques- 
tão central no estudo dos meios contínuos “mate- 
riais”. Conhecido este tensor, as equações de conser- 
vação: 


1) Ag TaB=0 


permitem, em princípio, determinar todo o evoluir 
dos corpos estudados (e do próprio espaço-tempo no 
caso do problema estar a ser abordado no âmbito da 
R.G., havendo, neste caso, naturalmente, que utilizar 
as equações de Einstein que incluem 1). 

Em ambos os livros indicados, os autores procu- 
ram estabelecer estas relações constitutivas estudan- 
do, principalmente Bressan, diferentes casos que in- 
cluem a transmissão do calor. Situam. muito jus- 
tamente, o problema no âmbito da Termodinâmica, e 
tiram do Segundo Princípio uma série de relações que 
condicionam as relações constitutivas admissíveis. 

Com surpresa nossa, porém, não encontramos nes- 
tes livros, em que há exaustivas referências bibliográ- 
ficas, qualquer referência às leis elásticas dos corpos 
rígidos válidas no caso adiabático, descobertas por 
processos extremamente simples pela primeira vez 
por Mc Crea (1952) [4,5], apresentadas num ante- 
rior artigo da Técnica [6] e que, sem dúvida, são uma 
primeiro exemplo preciso de relações constitutivas re- 
lativistas (3). 

São equações como as equações de Mc Crea que, 
em conjunto com as equações de conservação 1, per- 
mitem escrever as equações do movimento de um 
corpo sólido, como é o caso, por exemplo, da equa- 
ção da vibração de uma barra rígida, não homogénea 
apresentada em [6]. 

Não há nos livros citados sombra de uma equação 
do mesmo tipo, porque nunca é apresentado nenhum 
exemplo preciso de relação constitutiva que traduza a 
elasticidade dos sólidos (4). 

E se pode ser feito o comentário das equações de 
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Mc Crea dizerem respeito a um caso limite teórico, 
sem possibilidades de correspondência real, as fórmu- 
las elásticas dos sólidos não rígidos, por nós apresen- 
tadas em [7]. são bem as relações constitutivas rela- 
tivistas correspondentes à vulgaríssima lei de Hooke 
da Elasticidade Clássica. 

Vamos, nesta nota, usar estas fórmulas para esta- 
belecer as equações do movimento de uma barra não 
rígida. A concordância, na passagem ao limite c— 0, 
da equação obtida com as equações da Física Clássi- 
ca, atesta, a nosso ver, em favor da solidez e coerência 
das fórmulas e do método usado. 


Sendo (x', t) as coordenadas euclideanas e tempo- 
ral de um referencial S, o movimento de uma barra 
homogênea com uma lei elástica linear (lei de Hooke) 
verifica em Física Clássica a equação: 


2) a SR DT e 


dx? V2 dt? 


ra 


(x=x!,e X coordenada de Lagrange “fixa” à barra), 
invariante numa transformação de Galileu. 

Em Relatividade, o movimento de uma barra ho- 
mogénea rigida (não no sentido de indeformável mas 
no de mais rigida possivel) verificada a equação: 


d2x | DX 
à x? c2 dt? 


invariante numa transformação de Lorentz. (Ver, 
por exemplo [6]). 

A passagem de uma equação à outra não se faz 
pela simples passagem ao limite V— cc, porque a deri- 


2 w 2% 
vada — dá uma “volta” para - e passamos a ter 
BX< Ú Xx< 
02% : 32 x 
- onde anteriormente tinhamos 
o t- o t2 


Embora encoberta, a concordância entre a equação 
relativista e as equações clássicas é, no entanto, 
perfeita. 

Ão corpo rigido relativista que verifica 3 corres- 
ponde em F.C. o corpo rigido-indeformável que veri- 
fica 2 com V =00, 

Para estabelecer a relação entre a equação relativis- 
ta e as equações clássicas, há, pois, que escrever em 
primeiro lugar as equações relativistas dos corpos 
não rígidos. 

Utilizemos as fórmulas apresentadas em [7]: 


4) po=pº NI -F2+F2s2 -F2s 
Ri im 
(1 -F2)s 
É 


5) p=Po=o PF po -— 


| -F? 1 -F2+F2s2 


(po ep aa densidade da barra, deformada e não de- 
formada, no referencial próprio; s, relação entre os 
comprimentos de um elemento da barra, deformado e 
quando não deformado; F, coeficiente característico 
do material, O<F < 1) que, no caso limite de F= 1, 
nos dão as fórmulas de Mc Crea da elasticidade rela- 
tivista (a uma dimensão) dos corpos rigidos: 


o 
6) Pá = Po: gia <A 
Gê 
Cc 
Da 
7) P= Po- Po Edo =) 


Numa barra com o movimento em S dado por 
X = X(x, t), ou pela equação inversa x = x (X, t), 
temos: 


ox 
Te 
8) pe ME ac E 
t Y x 
X 
e 
Ô x 
| gesso via É 


| B? ias ] is é 
E e 7) 


As componentes T “ºBo do tensor impulsão-energia 
no referencial próprio local S, (que se desloca em 
relação a S com a velocidade v dada por 8) são: 


10) Taogo = 


o) 


 oõoS 
ooD o 
- 


P 
0 
O 
0 


Passando para as componentes do mesmo tensor em 
S obtemos: 


Po + B2 p,c? PotPoc? 
np Tél a ra 
0 O 0 0 
0 0 0 0 
Pot Be? 0 0 Bs cê Bº Po 
|- B? |-Bº 


Das equações de conservação | interessam-nos: 
12) daTId=0 e daT=O 


que nos dão, respectivamente: 


+BPpel). PotP,c2 o 
os 1-8? Cd E, 8) 
13) 
+ a D. Po + Pc? RM? sd 
+ 2 P) = ( = 


Feitas as substituições 4,5 e 9 e desenvolvidos os 
cálculos (tendo em conta que 8 é variável por o refe- 
rencial S, ser diferente de ponto para ponto e de 
instante para instante) as duas equações conduzem- 
nos à equação do movimento: 


2x | d2% da 
0x2 c2 dt? c2 c2F2 


= àx/ Ot2 
» OE VT dir 
bx dtax dtl. 


Chegamos assim por duas vias, a partir da conser- 
vação da quantidade de movimento (equação 1) e a 
partir da conservação da energia (equação 11), à 
equação do movimento de uma barra em coordena- 
das de Euler. 

Fazendo em 14 F =| 
do) obtemos 3. 

Para fazer a comparação com as equações da Fisi- 
ca Clássica convem-nos passar para uma representa- 
ção em coordenadas de Lagrange: x =x (x, t). 

O cálculo diferencial fornece-nos as fórmulas: 


(caso limite do corpo rigi- 


2. | 2x 
0x2 (3) 3 X2 
dx 

02% 0 d2x pEtqlejo= 02x 
. QRO! dxot XOX dx/ dt dXZ 
15) 

02X. DX dx (DX) /0XY Dx 

ot? dx dt? vã) (277) à x? 

dx XxX)? dx 
e (me) Te. 


que introduzidas em 14 nos dão a equação equivalente: 


— O o 


dx dx dx estã = |=0 


ox ot oxdt lox/ot? DL 


Fazendo a passagem ao limite c — cc e recordando as 
fórmulas: 


| | d 
17) F2 pº do (22) , 
“Ma ds |s=] Rs 


(ver [7]) obtemos a equação clássica 2. 


No caso de uma barra não homogénea F = F (X) e 


o 
= á 
ia Po (X), as duas vias, a da conservação da quan- 
tidade de movimento e a da conservação da energia, 


conduzem-nos à equação: 


dr do, 
18) Ae +B dx + C =0 


em que Céo primeiro membro da equação |4e A e B 
são as expressões: 


R-. 
ê 2(F? : ) 
Ea ad) === 
F2. | s2 F(F2 1 +) 
E peões 
(v F2 e pa =] |) 
s2 
e E 
Riad TE 
18”) B=- 2 


Fazendo F (X) = |, reencontramos, como caso par- 
ticular, a equação da vibração de uma barra rígida 
não homogénea apresentada em [6). 
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(1) Designamos neste texto por meios continuos “materiais” os 
meios em que é possível definir referenciais próprios locais. Fica 
excluído o campo electromagnético no vácuo. 

(2) No livro de Bressan encontramos indicação do uso da elastici- 
dade relativista no estudo dos instrumentos de detecção de even- 
ruais ondas de gravitação, assunto à que nos referimos em [7]. 


(3) Lichnerovicz no texto [8] obtem a equação de estado p = pº A 
[6] 
£ 


pic” para os líquidos “rígidos” (em que as ondas se propagam com 
a velocidade limite c), equação que decorre imediatamente das 
formulas de Me Crea também aplicáveis aos liquidos. O trabalhar 
com liquidos e não com sólidos impediu-o de chegar às fórmulas 
de base. 

(4) É feita referência a autores que, na falta de outros critérios, 
admitiram em Relatividade uma lei de Hooke linear, o que é des- 
mentido pelas fórmulas não lineares que apresentamos, 
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MODEL VALIDATION IN THE ESTIMATION OF 
NON-STATIONARY SPATIAL PHENOMENA 


HENRIQUE GARCIA PEREIRA, * Prof. IST 


ABSTRACT 


Experimental results of the application of a cova- 
riance (or variogram) model validation method for 
non-stationary phenomena are presented for the case 
of estimation of S grades in a pyrite deposit. 


RESUMO 


Apresentam-se os resultados experimentais da aplica- 
ção, à estimação do teor em S numa jazida de pirites, 
de um método de validação de um modelo de cova- 
riância (ou variograma) para fenómenos não-estaciond- 
rios. 


INTRODUCTION 


Many spatial phenomena ocurring in Nature may be 
described by means of Random Functions. Characte- 
ristic variables of these phenomena (measured in so- 
me experimental points) are viewed as the realization 
of a Random Function whose covariance must be 
inferred to allow the BLUE (best linear unbiased esti- 
mation) of the variable in unknown points or surfa- 
ces. The inference of that covariance from the expe- 
rimental data is not possible from one single realization 
of the phenomenon, when a stationarity hypothesis 
does not hold. 


For a non-stationary phenomenon, the traditional 
approach to deal with estimation problems is to split 
the variable Z(x) (where x are the coordinates of a 
point in one or two dimensions) in two terms: 


Z(x) =m(x) + Y(x) 


Where m(x) =E [Z(x) Jis the drift, an ordinary fune- 
tion that accounts for the regional component of the 
variable and Y(x) is the flutuation, a stationary ran- 
dom functiom with zero mean and stationary cova- 
riance that models the variability of the local compo- 
nent of the variable. 


If the covariance C(h) = E [FT Ex +h) Y(x)] was 
known, the problem of estimation of the variable in 
an arbitrary point x, could be solved on the grounds 
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of the Universal Kriging theory (MATHERON, 1970), 
using the system: 


AP C(Xa—X) =C(xa—Xo) +tHeft(Xo) a =1,2,..n 


[1] 


A“ Te(xa)=Tº(Xg) 6= 0, À É 


Where A“ and AÊ are the weights to apply to the 
values of the the variable in the n experimental points 
a in order to built up the estimator Z(Xo) ENO Lika Bio 
are Lagrange parameters and f* are monomials, de- 
pending on the degree k of the drift (m(x)= ac fe(x)). 
The Einstein summatión convention is used in the 


sequel. 


The kriging variance (minimum variance of Z(xo) — 
Z(x0)) is calculated by: 


o 4 = CHAO) tuctiXo) -A" Elis — Ko) [2] 


MODEL VALIDATION 


In order to apply system [1], a covariance model 
must be found. In Geostatistics (MATHERON, 1970), 
instead of covariance, it is worth to use o variogram 
function, as many phenomena display an infinite dis- 
persion. The variogram d (h) is calculated by: 


p(h)=VAR[Y(x+h)-Y(9)] 
=ElY(x+h)-Y (9) 


and system [1] holds if C(h) is substituted by — P(h). 


But the variogram function to which one has experi- 
mental acessis yp(h)= VAR[Z(x+h)-Z(x)] 
(variogram of the original variable Z(x)) and, except 
in some particular cases (PEREIRA, 1983), the va- 
riogram of Y(x) can not by calculated from data, 
because considerable biases are introduced. 


So, an a priori model for the variogram function 
must be set, and parameters of the model must be 
calculated. In this note, a method of model validation 


FIG. 1 — Experimental configuration of S grades (%) 


is applied to experimental data of S grades in a pyrite 
deposit (CVRM, 1978). 


The method consists in the computation of estimators 
Z(Xa ) (using Z(Xq )J=A8 Z(xg), YB%a) in all expe- 
rimental points where a value of Z(x q) exists. Those 
estimators are obtained from system [1], using the a 
priori model for y(h). The parameters of the model 
are choosen in a practical range whose limits are sug- 
gested by a preliminary analysis of experimental data 
(variograms of Z(x) and of Z(x) — m(x) beeing the 
basis of this analysis) and by all available qualitative 
information about the spatial development of the 
phenomenon. 


An utility surface is then built in function of two 
statistics: 


e The Average error À 
1 » * 
=| 2 (26) -Z0] 


e The Average Relative 
Quadratic error 


p= 2 [Z(X) = Z(009)]? Fo 


The shape of indiference curves in the neighbourhood 
of the optimum (given bye = 0,v = 1) depends on the 
relative “preference” for each kind of error (the first 
one,e, accounts for the average under or over estima- 
tion of the variable and the second one, v, measures 
the goodness of fit of kriging variances calculated 
from the model). 


The model choosen for further estimation procedures 
is the one for which values of statistics «e and v are 
closest to the optimum, in the utility surface space. 


EXPERIMENTAL RESULTS 


The above described method was applied to an expe- 
rimental configuration of 11 samples of pyrite ore 
assayed for S grades in a level of a pyrite deposit 
(Fis. 1). 

The variogram model tested was / (h) = wh,forwin 
the range [2,3 Jand fork = 0, 1,2 (degree of the drift). 


In the utility surface shown in Fig. 2 (each indifferen- 
ce curve beeing labeled by a code 1,4,6,9 denoting a 
growing deviation from the optimum e=0,v = 1), the 
experimental points for each pair of the model's pa- 
rameters w, k were plotted in function ofe and v. Itis 
worth noting that eis constant for each k value (pro- 
priety of linear variogram estimation) and values of «w 
are figures appearing in each straight line e =constant. 


FIG.2 — Plot of experimental points 
(«, k)in the utility surface space 
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An analysis of Fig. 2 allows the conclusion that va- 
lues «w = 2.6, k = O lead to a minimum deviation from 
the optimum (utility 4). 


Accordinglv, the variogram model to be used for further 
estimation (in the area depicted in Fig. 1) is y(h) = 
2.6 h with constant drift (k = 0). 
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A VARIABLE METRIC PROJECTION METHOD 
FOR MINIMIZING GENERAL FUNCTIONS 
SUBJECT TO UPPER AND LOWER BOUNDS 


ON THE VARIABLES 


MÁRIO J. A. LANÇA 


Centro de Electrónica Aplicada (INIC) e Dept. Eng. Elect. (Instituto Superior Técnico) 


RESUMO: 


Descreve-se um novo algoritmo de métrica varidvel 
para a minimização de funções de n varidveis sujeitas 
a limites. O algoritmo, que se destina a ser utilizado 
em programas para optimização de circuitos electró- 
nicos não-lineares, emprega a técnica do conjunto de 
restrições activas e tira partido da forma particular 
das restrições sobre as varidveis para simplificar os 
operadores matriciais envolvidos. Testes de compa- 
ração com outros métodos conhecidos forneceram 
indicações prometedoras. 


ABSTRACT: 


A new variable metric algorithm for minimization of 
general functions subject to upper and lower bounds 
on the variables is presented. The motivation came 
from the optimization of nonlinear electronic circuits 
to time-domain specifications. The constrained mi- 
nimization problem is treated through an active set 
strategy and full advantage of the particular form of 
the constraints is taken to simplify the manipulation 
of the various matrix operators. Comparative testing 
with published minimization procedures has been 
carried out with encouraging results. 


INTRODUCTION: 


In this note we consider the following nonlinear pro- 
gramming problem: 


Minimize f(x) (1) 
Subject to |; < x; < ujforj =1,2,...,n (2) 


where f(x) is a real-valued differentiable function of n 
continuous variables x. 


The motivation comes from the design optimiza- 
tion of nonlinear electrical networks to time-domain 
specifications. In fact, most design procedures are ba- 
sed on the minimization of a scalar performance in- 
dex which reflects a set of design specifications on the 
time-domain response. This index, the objective or 
error function, is typically an integral of the form 


E 
f(x) - | [Vo (t) -Y (OP dt, (3) 
Õ 
where V“(t) represents some desired output voltage. 
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The actual output response Vs (t) is evaluated by 
solving a system of nonlinear ordinary differential 
and algebraic equations which describe the circuit 
behaviour for the actual set of parameters X. 

It should be noted that box constraints on the de- 
signable parameters are invariably present in practi- 
cal problems and it is therefore desirable to handle 
these in an efficient way. The variable metric algo- 
rithm described in this note was specifically designed 
for this class of problems. 


OUTLINE OF THE MINIMIZATION PROCE- 
DURE: The specific aim of the algorithm is to locate 
a Kuhn-Tucker point, defined as a point at which f(x) 
is non-decreasing along any feasible direction. The 
search for this local constrained minimum is carried 
out through an active set strategy [1]: a certain num- 
ber k of constraints is assumed to be active at each itera- 
tion and the algorithm proceeds as 1Í the problem were 
one of unconstrained minimization in the subspace 
that corresponds to the intersection of the hyperpla- 
nes indexed by the active set. An exchange-type algo- 
rithm allows the active set to be varied from one itera- 
tion to the next. 

An iteration of the method begins with a test to 
assess the adequacy of removing a constraint from 
the basis. Here, we follow an approach due to Rosen 
[2] which is generally recognized as avoiding the zig- 
zagging phenomenon: a constraint should be remo- 
ved from the active set only if the gain thus likely to 
be achieved, with regard to a decrease in the function 
value, is estimated to be at least of the same magnitu- 
de as the reduction otherwise predicted. Furthermo- 
re, only those constraints with positive first - and 
second - order estimates of the Lagrange multipliers 
are considered as candidates to be dropped from the 
basis in order to ensure that both the steepest-descent 
and the Newton directions point into the feasible 
region. 

On each iteration, a local quadratic approximation 
(in the subset of the n-k currently active variables) is 
assumed which agrees with f(x) in value and gradient 
at the current point x; and whose Hessian is prescri- 
bed by a given symmetric matrix G|;: 


8: 6168 , (4) 


where 9=x —x; and gi; designates the projection of 


the gradient on the intersectior. =í the active cons- 
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traint hyperplanes. For a general function, 61; 1s an 
approximation to the projected inverse Hessian which 
is updated at each iteration according to the changes 
in x and g. The approximation (4) Is deemed to repre- 
sent satisfactorily f(x) over a region (hyperparallele- 
pipedon) 


Mx -xdl S hlu;- 1) (5) 


where h; is an appropriate scalar whose value is revi- 
sed at each iteration. 

Unlike Fletcher's method of hypercubes [3], our 
algorithm does not attempt to solve the constrained 
quadratic subproblem completely. Instead of this, a 
displacement of the form ô; =s; + 9 (v; — s;) is calcula- 
ted, where s; represents the distance to the constrai- 
ned minimum along the steepest-descent direction, 
Vi= G1 gi; designates the quasi-Newton correc- 
tion and O is the largest number in the range [0,1] 
that vields a feasible ô obeying (5). The directions s; 
and v; point into the feasible region because the algo- 
rithm ensures that the projected Hessian estimate re- 
mains positive definite and the rule for dropping 
constraints from the active set is such that first and 
second-order estimates of the Lagrange multipliers are 
positive. This strategy for selecting the displacement 
ô; 1s intended to bias the search towards the steepest- 
-descent direction during the first stages of the met- 
hod when é; (6) is not a close representation of f(x). 

The following rank-two correction formula due to 
Powell [4] is used to update the Hessian estimate 


ui na 
G'-G+a 48 +ô8u q2dô (u!0) (6) 
[8º “ot 


where yu =y — Gô and y-g(x;j+ ô;) -g(x;). No 
more than one negative cigenvalue can be introduced 
by this update for a e [0,1]. It is therefore suffi- 
cient to choose a damping factor « in this interval 
vielding det G'> O to ensure that the Hessian estimate 
remains positive definite [5]. 

The algorithm also keeps and updates at each itera- 
tion estimates of the inverse Hessian, D, and of the 
projected inverse Hessian in the (n-k) — dimensional 
subspace of the variables that are not indexed by the 
current constraint basis, Hj = G4. This operator is 
modified by means of appropriate recurrence rela- 
tions each time the active set is altered. These rela- 
tions are such that the information gained in previous 
iterations about the inverse curvature of f(x) along 
the direction orthogonal to an hyperplane dropped 
from the basis is contained in the new operator H,. 

The updating of the inverse Hessian estimate D is 
performed by means of the following correction 
formula: 


D'-D [ue Db) -(n5! D+ Dan”) (87 Dy") +(D66" DIn!y 1] (7) 


Li 
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whereo =(uTy) (6TD6) -(8TDy), withD=G"' 

andmn = Du=Dy 8. The damping effect of a is 

simulated in (7) by means of the modified vector 
a(l-a) (6! u) 


ÉS Djs O & (8) 


The projected inverse Hessian estimate H | is updated 
by the (n-k) — dimensional equivalent of (7). 


ALGORITHM PERFORMANCE: The testing of 
the algorithm has been done against three well known 
unconstrained minimization procedures: the Davidon- 
-Fletcher-Powell algorithm [6], the Powell's trajec- 
tories method [4 ]and the Fletcher-Gill-Murray met- 
hod [7]. The transformation of variables y; =1; + 
(u;-1;) sin? x; was used to accommodate the box 
constraints when employing these procedures to solve 
constrained problems. The results of the test on three 
test functions (a quadratic function, Rosenbrock's 
parabolic valley [8] and Wood's function of four varia- 
bles [9]) provided the indication that the new algorithm 
performs consistently better than the other three in 
those problems which present a minimal point lying 
on the boundaries, and not worse, on the average, for 
the other cases. 

Practical circuit design problems have also been 
satisfactorily solved by the algorithm [5 ]. 


CONCLUSIONS:To sum up then, a variable metric 
algorithm for minimizing a function of several varia- 
bles was described. The algorithm was tailored to suit 
the needs of automated design of electronic circuits. 
Numerical evidence was referred showing the effi- 
ciency of the algorithm in the number of function 
evaluations required to solve well known test problems. 
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THE ESTIMATION OF TRANSFER FUNCTIONS 
IN NATURAL RESOURCES EVALUATION 


FERNANDO DE OLIVEIRA MUGE, Prof. Auxiliar do IST 


ABSTRACT: 


Experimental results concerning the estimation of 
direct and indirect transfer functions are presented 
for the case of an iron reboly. 


RESUMO: 


São apresentados resultados experimentais da es- 
timação de funções de recuperação directas e indirec- 
tas para o caso de um jazigo de ferro. 


INTRODUCTION: 


One of the problems that arises in the exploitation 
of a natural resource is its quantification in terms of 
some selection or cut-off parameters. In what con- 
cerns, for instance, the exploitation stage of a mineral 
resource, it is very important to forecast the grade- 
-tonnage curves which are a basic tool in the study of 
the feasibility of a mining venture. These curves are 
generally called Transfer Functions (MARECHAL, 
1974; MATHERON, 1974). 

In the estimation of these Transfer Functions (T.F.) 
two features must be taken into account: 


e The support of the selection is a spatial dominion 
v(block). 


e Every selection criterion will be applied to estima- 
” * - à 
tiors Z of the average value of the regionalised 


variable (R.V.)Z(x) in the dominion v, given the 
available information (Z(x,). «q = |,...n)inthen 
experimental points. 


The Direct Transfer Functions (D.T.F.) take into 
account the support effect and the Indirect Transfer 
Functions (I.T.F.) the so-called information effect 
(MUGE, 1979). 


THE DISCRETE GAUSSIAN MODEL 
The essential assumptions of this model are: 


[D Asa starting point the global field in which it 
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Is spread the natural resource is conceived as 
the union of a finite number of small blocks v. 
A random value Z(x) s associated with each v 
so that it gives a discrete version of the moving 
average of the punctual stationary Random 
Function (R.F.) Z(x). 


[) The R.F. Z(x) is considered as a Gaussian 
Transform ó (Y(x)) of a Gaussian R.F. Y(x) 
and Z, as the Gaussian Transform 4, (Y,), where 
Y,is also a Gaussian R,F. 


HI) It is assumed a permanece of law hypothesis 
which is used to link the parameters of the two 
Gaussian Transform: gand d.. 


IV) The two dimensional law of the pair (Y(x), Y,) 
is bivariate normal with a correlation coeffici- 
ent r > o, independent of x. 


In the frame of this model the expression of dy (Y,) 
becomes: 


> Cn 
dv(Yy) = n 20 nt PH (14) (1) 


with Cn= f & (x) Hn(x) g(x) dx being the Fourier 
coefficients of the expansion of & in Hermite Polyno- 
mials, 


The estimation of a D.T.F. like the tonnage-cut-off 
grade curve (T(Z,)) becomes: 


TZ) =1 =6(Y;) (2) 


Where: Z, = &v, (Y,) and G(Y,) is the cumulative 
gaussian distribution function. 


If the estimator Z. of Z, is of the form: 
Z =% ha Zig) with & hos | (3) 
where Ay are the weights given to the values of the 
R.V. in the n experimental points, it is possible to 
apply the same model for the estimation of 1. T.F.. 
EXPERIMENTAL RESULTS 

The above described model was used to estimate 


some T.F. of an iron orebody. 
Firstly the model was checked by comparing the 
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FIGURA 1 — Experimental and Theoretical Histograms 
of Sample Grades 


experimental histogram of the 791 core samples that 
covers all the deposit with the theoretical one dedu- 
ced from the adjusted gaussian transform (Fig. 1). 

As it is observed, not only the mean and variance 
are very close to the experimental ones, but also the 
theoretical histogram fits quite well. 

Another check consists in the estimation of some 
T.F. based upon sample values (nearly point support). 
In Fig. 2 are presented the point support theoretical 
and experimental tonnage-cut-off grade curves. The 
theoretical curve matches perfectly the experimental 
one. 


FIGURA 2 — Point Support Theoretical and Experimental 
Tonnage — Cut — Off Grade Curves 


Once adjusted and checked, the model was used to 
estimate the histogram and the tonnage-cut-off grade 
curve (I.T.F.) of a future selection procedure based 
upon 50x25x10m block estimates. 

The estimator envisaged for these blocks was a 
BLUE (best linear unbiased estimator) or kriging es- 
timator (MATHERON, 1970). 

The resulting histogram is very different from the 
experimental one based sample values: the varian- 
ce is almost 7 times less and it is much more smooth 
(Fig. 3). 

The effect of using the point support tonnage-cut- 


FIGURA 4 — Point Support Experimental Tonnage — Cut 
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FIGURA 3 — Experimental Histogram of Sample Grades 
and Theoretical Histogram of 50x25xIOm 
Block Estimates 


-off grade curve (sample values) instead of the corres- 
ponding 1.T.F. will be the over estimation of the ton- 
nage in the zones of low and high cut-off grade (Fig. 4). 
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— Off Grade Curve and I.T.F. of 50x25x10m 
Block Estimates 
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anuncie 
na técnica 


Com uma tiragem de 5 000 exemplares e uma difusão 

à escala nacional e internacional a “TÉCNICA” abarca nos seus leitores 
quadros científicos e técnicos das mais variadas empresas, bem como 
grande número de professores e estudantes universitários. 

Face a isto, poderemos concluir que se a “TÉCNICA” 

merece o apoio de todos os que de algum modo se encontram 
ligados à Engenharia — Escolas, Institutos de Investigação e Empresas 
— ela é, para além de um veículo da investigação teórica e experimental, 
um veículo por excelência de produtos e serviços, em suma, 

das actividas da vossa empresa. 

Através da inserção de publicidade da vossa empresa na revista 
“TÉCNICA” está apoiando esta publicação e comunicando 

com um universo de potenciais interessados nas vossas actividades. 
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